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ANNALEN DER PHYSIK 


FOLGE * BAND 17, HEFT 1 x 1955 


Zur wellenmechanischen Behandlung 
von Vielkörperproblemen 


Von Wilhelm Macke 
Mit 2 Abbildungen 


Inhaltsübersicht 


Ein Verfahren zur quantenmechanischen Behandlung vieler Fermiteil- 
chen wird behandelt, dessen nullte Näherung die statistische Theorie und ihre 
Verfeinerungen liefert. In höherer Näherung approximiert dieses Verfahren 
die Hartree-Focksche Methode zur Berechnung solcher Vielkörperprobleme. 


Einleitung 


Physikalische Probleme, an denen viele Teilchen beteiligt sind, die der 
Fermistatistik genügen, lassen sich unter Vernachlässigung der Correlations- 
wirkungen nach der Methode von Hartree und Fock behandeln, bei der die 
wahre Wechselwirkung der Teilchen durch eine mittlere Wechselwirkung er- 
setzt wird und eine Beschreibung durch Teilchenfunktionen %, (t) möglich ist. 
Das Pauliprinzip bleibt durch die Forderung nach höchstens zweifacher Be- 
setzung der Zustände y, gewahrt, die untereinander ein orthogonales System 
von Funktionen bilden. 

Das Verfahren der Berechnung der y, nach der Methode des „self con- 
sistent field‘ ist sehr mühevoll und bei großer Teilchenzahl N praktisch un- 
durchführbar. Für eine genäherte Bestimmung der Wellenfunktionen y, er- 
gibt sich dagegen die Schwierigkeit, daß diese Näherungsfunktionen im allge- 
meinen kein Orthogonalsystem bilden, so daß die Erfüllung des Pauliprinzips 
fraglich ist und zum Beispiel Variationen der Energie nach irgendwelchen 
Parametern der y, nicht durchgeführt werden dürfen, da bei der Variation 
möglicherweise alle Teilchen in den Grundzustand ‚rutschen‘ könnten. Aus 
diesem Grunde hat man bislang bei Problemen mit großer Teilchenzahl auf 
eine wellenmechanische Behandlung verzichtet und statistische Beschrei- 
bungen aufgesucht, bei denen immerhin das Pauliprinzip gewahrt bleibt!) ?). 
Eine Reihe von Arbeiten befaßt sich darüber hinaus mit möglichen Verfei- 
nerungen des statistischen Modells*)*). Da die Herleitung dieser Verfeinerun- 
gen aber nicht über die Angabe von Wellenfunktionen erfolgt, bleibt ihre 
quantenmechanische Basis unsicher. 


') L. H. Thomas, Proc. Cambridge philos. Soc. 28, 542 (1927). 

*) E. Fermi, Z. Physik 48, 73 (1928). 

8) C. F. v. Weizsäcker, Z. Physik 96, 431 (1935). 

*) P. Gombäs, Acta Physica Hungaria 3, 105 (1953); 3, 127 (1953). 
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Aus den angeführten Gründen wird hier ein Näherungsansatz der 
Wellenfunktionen y, für große Teilchenzahlen N behandelt, bei dem die 
y,„ Stets ein Orthogonalsystem bilden. Alle Funktionen y, enthalten dabei 
noch eine oder mehrere willkürliche Funktionen, die durch Variation der 
Energie bestimmt werden. So entsteht ein Näherungsverfahren, das in nullter 
Näherung die Gleichungen des statistischen Atoms liefert, das sich unter ent- 
sprechendem Rechenaufwand beliebig verbessern läßt und durch das die 
Hartree-Focksche Methode sukzessiv approximiert werden kann. In 
dieser Arbeit ist der Schwerpunkt auf die Reproduktion der statistischen 
Theorie gelegt, während die durchfihrbaren Verbesserungen erst am Schluß 
diskutiert werden und den Inhalt einer weiteren Arbeit bilden. 


Dichte und kinetische Energie 

Beim exakten N-Körperproblem ist die Energie E ein Funktional der 
Wellenfunktion Y, und ihr Wert bestimmt sich durch Variation der Energie 
über alle möglichen, aber normierten Wellenfunktionen: 


Die Näherung statistisch unabhängiger Teilchen geht aus diesem Ver- 
fahren (1) hervor, wenn man für die Wellenfunktion Y eine Determinante aus 
N Teilchenfunktionen y, (k) ansetzt 


W1...N)= 


1 
In diesem Falle hängt E von den Funktionen y, einzeln ab und (1) geht über 
in das Verfahren von Hartree und Fock: 


OE ¢ 
E == E [yı Wo .. oy, 0 (Pp Ym) Ö, m’ (3) 


(k) |. (2) 


Teilchendichte 9 und kinetische Energie haben in dieser Näherung die Form 


o= > Exin = Jat > : (4) 


n=1 n= 


soll durch spezielle Voraussetzungen über die y, genähert gelöst werden. 

Der eigentlichen Behandlung voraus geht eine formale Umformung der 
Ausdrücke (4) für o und Ex;„, die noch keinerlei Einschränkung für die Funk- 
tionen y, bedeutet. Es wird eine zunächst willkürliche, reelle Funktion oft) 
eingeführt durch den Ansatz 


Ym (t) = Pm (5) 


für die y,. Diese Gleichung kann bei gegebenem o zunächst als Defini- 
tionsgleichung von @,, angesehen werden. Die Teilchendichte o wird 


N 
(6) 


also eine Funktion von o und den 9, 


Das Problem (3) der Berechnung von Energie E und Teilchenfunktionen y, | 
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Für die kinetische Energie Ex; aus (4) läßt sich ein Ausdruck in o, o und 
den p, angeben, der für das Weitere von besonderer Bedeutung ist. Für das 
Quadrat des Gradienten einer Funktion Vor mit reellen o und eventuell 
komplexem r läßt sich folgende Formel angeben: 

öyor öyr ? 1 da 

Diese Formel wird zunächst für t = o/o verwendet und u 
öye 3 Volo 3 1 da da/o 0 2 


or et 2a 


Die Anwendung der gleichen Formel (7) mit |/ot = y,, also Yr = 


1 
VN Pm 


ergibt in der Summe über die n: 


Die beiden letzten Terme von (8) und (9) sind miteinander identisch wegen (6). 
Mit den Gin. (6), (8) unf (9) kann deswegen die kinetische Energie (4) in der 
Form 

Exin = farts "WW +( or ox 
dargestellt werden. Diese 3 a der kinetischen Energie werden später die 
Fermienergie, den Weizsäckerschen Energieanteil und ein für große N 
verschwindendes Zusatzglied beschreiben. Vorläufig stellen (6) und (10) 
bloße Umformungen der Dichte o und der kinetischen Enegrie Ey; dar auf 
Grund des Ansatzes (5) für y,. 


Der Variationsansatz 


Für die Wellenfunktionen y, soll nun ein Ansatz gemacht werden, der 
eine genäherte Auswertung des Problems im Sinne der Gleichungen (3) er- 
möglicht. Zu diesem Zweck wird irgendein bekanntes, aber sonst beliebiges 
Orthogonalsystem von Funktionen ¢ (1) ‘acest die die Relationen 


Jat, (9) = (11) 
voraussetzungsgemäß erfüllen, wenn die ion über den gesamten 
Variationsbereich der Koordinaten 4 läuft, in dem dieses Funktionensystem 
erklärt ist. 


Zwischen den Raumpunkten r und den Punkten y sei eine eindeutige, aber 
sonst willkürliche Transformation gegeben: 


y = y(t) | (12) 


über die später in speziellererer Weise verfügt wird. Nun läßt sich für die 
Funktionen y, folgender Näherungsansatz machen: 
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Das durch (13) gegebene System von Funktionen yp) (r) ist ein Orthogonal- 
system, denn es ist 


(vn. ym) = dr Gn = J Gr. (Y)* Ym (9) = On m (14) 


entsprechend (11). Die Funktionen g!, sind bekannt. Der Ansatz (13) für alle 
N Funktionen y}, enthält also nur die drei willkürlichen Funktionen 4(t), 
über die später so verfügt werden kann, daß die Energie ein Minimum wird. 
Die Orthogonalität der y), bleibt bei dieser Variation gewahrt. Im übrigen hat 
der Ansatz (13) genau die Form (5), so daß Dichte o und kinetische Energie E, 
durch (6) und (10) dargestellt sind. 


Die statistische Näherung 


Zur Auswertung der kinetischen Energie müssen spezielle Orthogonal- 
systeme 9, und Transformationen herangezogen werden. Will man lediglich 
die Energie (3) in der statistischen Näherung N—0oo bestimmen, so ist die 
spezielle Wahl der Funktionen 9), nicht sehr wichtig und kann durch Zweck- 
mäßigkeitsgesichtspunkte bestimmt sein. Bei Verfeinerungen der Theorie 
werden allerdings die Ausgangsfunktionen 9), eine wichtigere Rolle spielen, 
und man wird solchen Funktionen den Vorzug geben, die der Kugelsymmetrie 
des Problems Rechnung tragen. 

Hier werden zunächst unter den Transformationen 4) (r) nur solche zuge- 
lassen, die den Raum r in einen Einheitswürfel transformieren 


sy, & = 1, 2,3 (15) 
und deren Flächenscharen y, — const ein Orthogonalsystem bilden: 


OY; OYR 
er ör Dix hi, (r). 


Dann existiert nämlich für a die Darstellung 


eo=N 


= NViheh- (17) 


OX, 


Cx, or, 


Als Funktionen 5 werden die Eigenfunktionen des Einheitswürfels 
gewählt, also 


0 


3 
gn = II 2sin yp (18) 


Unter diesen Voraussetzungen lassen sich (6) und (10) für große N weitgehend 
vereinfachen. 

Die Funktionen (18) haben die Eigenschaft, bei gleichmäßiger Besetzung 
aller unteren N Zustände, für große N die Relation 


1 


N->oo N (19) 
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zu erfüllen. Man macht sich diesen Sachverhalt (19) leicht am eindimensionalen 
Beispiel klar, das in Abb. 1 dargestellt ist. Im Mittel, das heißt bis auf Schwan- 
kungen, die mit wachsendem N klein F 

werden, konvergiert die Summe gegen 
N, der Gesamtausdruck (19)also gegen 1. 
Aus (6) und (19) folgt unmittelbar: 15 


0>o. (20) 


o bedeutet also, abgesehen von örtlichen 
Schwankungen, die mittlere Dichte o 
und y(t) in (12) transformiert die 95 
Koordinaten r gerade so, daß in den 
neuen Koordinaten die Dichte im Mittel 
konstant, nämlich gleich 1 ist. 
In der durch (18), (19) und (20) ge- 0 

gebenen Näherung großer Teilchen- Abb.1. Darstellung der Funktion f(y) 
zahlen ist es möglich, die kinetische wit: rz (V2 sinany) fir N =7. Die 
Energie (10) als Funktion der Dichte o 
bzw. o allein darzustellen. Der dritte 
Term in (10) verschwindet unter der 
Voraussetzung (20), also o > o, der zweite hängt ohnehin nur von gab, während 
der erste, im weiteren mit E7 bezeichnet, noch einiger Umformungen bedarf. 


1.0 


ei a von 1 werden mit zu- 
nehmendem N kleiner 


Auswertung der kinetischen Energie 
In diesem Abschnitt soll der erste Term der kinetischen Energie (10) 


im Sinne des Ansatzes (18) fir (pn ausgewertet werden. Aus der indirekten 
Abhängigkeit der g?, von r folgt zunächst 


Dieser Ausdruck vereinfacht sich wegen der vorausgesetzten Orthogonalität (16) 
der Transformation t) (tr), und weiter bei Bezugnahme auf die spezielle Form (18) 
der Funktionen g), zu: 


Jar > mi (Pn)? (23) 
n k 


Bei großen Teilchenzahlen fallen in (23) unter den (g))? hauptsächlich die 
Funktionen mit großen Indizes n,, ny, n, ins Gewicht. Diese Funktionen 
2 sin? n, y, oszillieren stark zwischen 0 und 2 in Bereichen, in denen o und 
f, als langsam veränderliche Funktionen praktisch konstant sind. In der 
Näherung N — oo können daher die (9)? durch ihre Mittelwerte, also durch 1, 
ersetzt werden. 
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Nun bleibt noch der Ausdruck 


[A (3) (3) (3) ] (24) 


Ny, Nay Ny 


auszuwerten. Mit Einführung von 


m=1,2,3,... (25) 


bekommt die Summe in (24) die einfache Form 
 pe+ 22 + 22] = ve. 26 
=[A(3) ] = [ 2 4 — ( ) 


Stellt man sich einen dreidimensionalen Raum mit den Koordinaten k,, k, und 
k, vor, so werden durch (25) in diesem Raum alle Punkte eines Gitters ausge- 
zeichnet, und zwar in demjenigen Achtelraum, in dem alle drei Koordinaten 


positiv sind. Die zugehörigen drei Gitterkonstanten sind N, fe und I 


Die Summation in (24) bzw. (26) erfolgt also über alle positiven Gitterpunkte 
des f-Raums, während das Abstandsquadrat f? jedes Punkts vom Koordinaten- 
ursprung einen Beitrag zur Energie 
liefert. 

Bislang ist noch nicht festgelegt, 
welche dieser Punkte als ‚besetzt‘, 
also als in der Summe enthalten, an- 
zusehen sind. Da voraussetzungsgemäß 
der niedrigste Energiezustand des Sy- 
stems aufzusuchen ist, werden in (24) 
und (26) diejenigen N Gitterpunkte zu 
jE besetzen sein, die zur Gesamtenergie 
$ Mr am wenigsten beitragen, deren f? also 

möglichst klein ist. Die im Grund- 
Abb. 2. Darstellung des t-Raumes in Zustand des Systems besetzten Gitter- 
zwei Dimensionen. Unabhängig von der punkte liegen daher alle in einer Achtel- 
verschiedenen Größe der Gitterabstände . kugel, siehe Abb. 2, deren Radius ge- 
liegen emer rade so groß ist, daß er die zentrums- 
nächsten N-Punkte umfaßt. 

Bei der mathematischen Durchführung der Überlegungen des letzten Ab- 
satzes können in der Näherung großer N die Summationen durch Integra- 
tionen ersetzt werden. Zunächst ist der Radius X der Kugel zu bestimmen. 
Bei g = 2facher Besetzung jedes Zustandes (g = 2s + 1, s = !/,) lautet diese 
Bedingung: 


vy = 1= + San, dn, dng = 1 (27) 


By Ng My 


und nach Ersetzung der n, durch die k, gemäß (25) und Integration über die 
Achtelkugel 


Vis 


K 
89 
=~, dk, dk, dk, = — | 4ak* dk = Ks, 28 

hi! 1 2 3 =F 30 ( ) 


Dabei ist Gebrauch von (17) gemacht. 
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Mit diesem Wert für K läßt sich die Summation (26) bestimmen: 


5 5 \4ag/ 


Die Fermische Energie Ey wird also mit (24), (26) und (29) sowie o > o 


Der Gesamtausdruck für die kinetische Energie wird so 


Exin fee Ar (eVe) | (31) 


Er enthält also gerade die Energieterme von Fermi und Weizsäcker und 
hängt überhaupt nicht mehr von der speziellen Transformation  (r) ab, sondern 
über o und o nur noch von der Funktionaldeterminanten dieser Transformation. 

Genau genommen, wäre in (21), (23) und (24) erst die Integration durch- 
zuführen und dann die Summation über die Zustände n. Dabei würde dann die 
Wurzel in (28) nicht f,, fg, fs, sondern entsprechende Mittelwerte f,, fa, /, ent- 
halten und die Energie daher etwas größer werden. Der Unterschied zwischen 
beiden Verfahren liegt nämlich darin, daß in (28) die Wahl der N zu besetzen- 
den Zustände schwach, nämlich mit den f,, ortsabhängig gewählt ist. Man über- 
zeugt sich leicht, daß das Pauliprinzip dabei stets gewahrt bleibt. Vielmehr 
entspricht das in (27) bis (31) durchgeführte Verfahren einem allgemeineren 
Ansatz als (13) für die Wellenfunktion, nämlich 


statt yh, (32) 


m 


wobei in der Auswertung über die schwache Ortsabhängigkeit der a,, so ver- 
fügt wurde, daß die kinetische Energie zum Minimum wurde. Im übrigen ist 
der Unterschied zwischen beiden Verfahren bei großen Teilchenzahlen gering, 
da nur eine äußere Schicht der Kugel von dieser zusätzlichen Variation be- 
troffen wird. 


Variation der Teilchendichte 


Zur Berechnung von E und o ist es erforderlich, die Gesamtenergie durch 
Variation der Wellenfunktionen bzw. deren Parameter zu einem Minimum zu 
machen. Es wird für das Folgende vorausgesetzt, daß die potentielle Energie E ,, 
genau wie vorher die kinetische Energie, durch ein von o allein abhängiges 
Funktional darstellbar sei. Als Variationsableitung von E, wird dann die 
Größe Ufo] eingeführt. 

OE 

E, = E, [el to =U fe]. (33) 
Bewegen sich zum Beispiel die Teilchen in einem gemeinsamen äußeren 
Potential V(r), so ist 


Ufo = Vi). (34) 
Für ein zwischen je zwei Partnern wirkendes Potential V(r — r’) gilt 
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Die zugehörige, von der gemischten Dichte abhängige Austauschenergie läßt 
sich für große N im allgemeinen ebenfalls als von o bzw. o allein abhängig dar- 
stellen. Für die elektrostatische Austauschwechselwirkung zwischen Elek- 
tronen gilt unter Verwendung bekannter Näherungsverfahren zum Beispiel 


2 ‘ 


(36) 


Selbst die an sich bei diesem Verfahren nicht berücksichtigte Korrelations- 
energie®)®) läßt sich in der Form (33) genähert berücksichtigen’). 
Die Gesamtenergie des betrachteten N-Teilchensystems ist somit 


aVa\2 
E= + (?te) | + E, [ol (37) 


und hängt in der betrachteten Näherung nur noch von o baw. 9 ab. Will man, 
um E und 9 zu berechnen, über die freien Parameter der durch (13) gegebenen 
Wellenfunktionen variieren, so ist dabei zu berücksichtigen, daß die Gesamt- 
energie E in (39) gar nicht mehr von der willkürlichen Transformation (12), 
also y (r), direkt abhängt, sondern nur noch von deren Funktionaldetermi- 
nante oa ~ o. Das Volumen des Bereichs (15) liefert für o die Bedingung 


fdco=N fac = N San N, (38) 


die mit der Normierungsbedingung fiir 9 identisch ist. 
Es ist daher nur noch über o —o zu variieren unter Wahrung der Nor- 
mierungsbedingung (38). Diese Variation hat die Form 


6 {E—V, f dro} = 0 (39) 
und führt unmittelbar auf die bekannte Bestimmungsgleichung 


i? 2 r r ER 
+ +U =0, (40) 


während V, sich aus der Normierungsbedingung (38) fiir 9 bestimmt. 


Anwendungen 


Die Gln. (37) und (40) sind somit auf wellenmechanischer Basis abgeleitet. 
und es bleibt die Frage zu diskutieren, in welcher Richtung sich die durchge- 
führten Überlegungen verbessern lassen. 

Die erste und wichtigste Verbesserung besteht darin, daß mit y (r) die 
Koordinaten r nicht auf einen Einheitskubus, sondern auf eine Einheitskugel 
transformiert werden. Für die Wellenfunktionen p, sind dann die Eigen- 
funktionen der Potentialkugel anzusetzen, bestehend aus Kugelfunktionen. 
Yim (9,9) und halbzahligen Besselfunktionen (R), wobei die gesamte 
Transformation durch R (r) gegeben ist. Die Funktion o steht mit R(r) in 
5) E. Wigner, Physic. Rev. 46, 1002 (1934). 

6) W. Macke, Z. Naturforsch. 5a, 192 (1950). 
7) H. Hellmann, Quantenchemie, Leipzig 1937. 
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einem einfachen, umkehrbaren Zusammenhang, und die Dichte @ in (6) zer- 
fällt — voll besetzte Schalen vorausgesetzt — in eine Summe 


e= (r) (41) 


von Gliedern o,, deren jedes die Verteilung der Teilchen mit dem Drehimpuls / 
für sich enthält, und die später unabhänig voneinander variiert werden können. 
An die Stelle von (40) tritt ein System von Bestimmungsgleichungen für die 
einzelnen @,, deren eventuelle Verkopplung miteinander von der speziellen 
Form der vorausgesetzten Wechselwirkung abhängt. 

Der nächste Schritt zur Verfeinerung der Rechnung besteht darin, die 
Voraussetzung großer N nur noch bedingt aufrechtzuerhalten und die nächsten 
Glieder einer entsprechenden Entwicklung für die Bestimmung von o und 9, 
mit zu berücksichtigen. Man erhält so genauere, allerdings auch entsprechend 
kompliziertere Bestimmungsgleichungen für o wie für die o;. 

Schließlich läßt sich das Verfahren auffassen als mullte Näherung eines 
Schrödingerschen Störungsverfahrens zur Bestimmung der exakten Har- 
tree-Fockschen Wellenfunktionen (5) 


Yn = Yat Yn (42) 
Unter wi in (42) ist dabei eine der hier verwandten Wellenfunktionen (13) zu 
verstehen. Die in (13) noch enthaltenen Funktionen o bzw. R(r) seien durch 
(40) bzw. eins der anschließend diskutierten Verfahren bestimmt. Durch Ver- 
gleich der dann für y) und der für y, gültigen Bestimmungsgleichungen läßt 
sich schließlich ein Ausdruck auffinden, der in dem Störungsverfahren (42) 
die Rolle einer Schrödingerschen Störungsenergie übernimmt. 

Die erste Näherung dieses Verfahrens (42) besteht darin, die Energie mit 
dem exakten Hamiltonoperator, aber den Wellenfunktionen y), nullter 
Ordnung zu bestimmen. In dieser ersten Näherung macht sich die Schalen- 
struktur der Fermigase bereits qualitativ bemerkbar, wie durch entspre- 
chende Abschätzungen gezeigt werden konnte. Unter sehr allgemeinen Voraus- 
setzungen ergeben sich hier die richtigen, durch die Schalen ausgezeichneten . 
Zahlen sowohl für die Atomhülle als auch entsprechend für den Kern, für 
letzteren allerdings unter der üblichen Voraussetzung der Existenz einer 
starken Kopplung zwischen Spin und Bahn der Nukleonen. 


Dresden, Institut für theoretische Physik der Technischen Hochschule. 


Bei der Redaktion eingegangen am 30. April 1955. 


Ann, Physik. 6. Folge, Bd. 17 


ar- | 
viel 
36) 
ns- 
37) 
an, 
en 
nt- 

2), 

ni- 

ing 

38) 

or- 

39) 

4()) 

et. 

ge- | 

die 

gel 

len. 

ate 

in 
| 1b 


Nichtstationäre Vorgänge in Photoleitern. I 


Langzeitiges Anklingen der Photoleitung 
in CdS-Einkristallen als Methode zur Störtermanalyse 


Von K. W. Böer und H. Vogel 
Mit 7 Abbildungen 


Inhaltsiibersicht 


Das langzeitige Anklingen der Photoleitung bei schwacher optischer An- 
regung wird genauer untersucht und an Hand allgemeiner reaktionskinetischer 
Modelle diskutiert. Es zeigt sich, daß das Auftreten von Wendepunkten in 
der Anklingkurve notwendig an die Existenz von Akzeptortermen gebunden 
ist und daß aus der Lage dieser Wendepunkte Aussagen über die Konzen- 
tration der Akzeptortermgruppen und die entsprechenden Übergangswahr- 
scheinlichkeiten einfach hergeleitet werden können. 


Einleitung 


Die Betrachtung nichtstationärer Vorgänge ist in der Festkörperphysik 
von großer Bedeutung. So können auf diesem Wege die energetische Verteilung 
von Störtermen, die Besetzung dieser Terme mit Elektronen sowie vor allem 
die Übergangswahrscheinlichkeiten für Elektronen aus diesen räumlich lokali- 
sierten Zuständen in angeregte lokalisierte (exeited states) oder nichtlokali- 
sierte Zustände (Bänder) näher untersucht werden. Gewöhnlich werden zu 
solchen Untersuchungen die Lumineszenz oder die elektronische Leitfähigkeit 
herangezogen ; ihr An- und Abklingen nach unstetiger Anderung der Anregung 
(z. B. Ein- und Ausschalten), Glowvorgänge u.ä. werden dabei näher 
analysiert. 

Zur Auswertung dieser Untersuchungen sind ,,reaktionskinetische Modelle“ 
zugrunde zu legen, in denen die zeitliche Anderung der Elektronenkonzen- 
tration in jeweils einer Gruppe von Zuständen (z. B. Störtermen, Bändern 
usw.) durch die Differenz aus Zu- und Abgängen beschrieben wird. Schon 
relativ einfache Modelle aus nur wenigen Gruppen solcher Zustände!) führen 
auf Systeme simultaner nichtlinearer Differentialgleichungen, die nicht ge 
schlossen lösbar sind. Zudem enthalten solche Systeme eine relativ große 
Anzahl im allgemeinen unbekannter Parameter, z. B. Stértermkonzentra- 


!) Z. B. auch ein in letzter Zeit etwas genauer untersuchtes System aus zwei Bändern 
he und Valenzband) und einer Störtermgruppe (Akzeptoren) (vgl. z. B. Adiro- 
witsch?)). 

*) E. I. Adirowitsch, Einige Fragen zur Theorie der Lumineszenz der Kristalle, 
Deutsche Ubersetzung: Akademieverlag Berlin 1953. 
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tionen und Übergangswahrscheinlichkeiten, so daß auch eine numerische 
Lösung nicht sehr sinnvoll erscheint. 

Es ist nämlich von vornherein keineswegs gesichert, daß bei der recht mühsamen 
numerischen Integration eines solchen Differentialgleichungssystems durch Anpassung 
an den experimentell gefundenen Verlauf eine eindeutige Bestimmung der unbekannten 
Parameter möglich ist. 

Umgekehrt kann man versuchen, die experimentellen Ergebnisse als Lö- 
sungkurven sehr einfacher Differentialgleichungen darzustellen. Für jeweils 
eine Differentialgleichung geht dies jedoch erfahrungsgemäß nur in begrenzten 
Bereichen der untersuchten Vorgänge (z. B. in gewissen Zeitbereichen). In 
dem Bereich, in dem die Lösungskurve der angenommenen Differentialglei- 
chung mit der experimentell gefundenen Kurve übereinstimmt, kann diese 
Gleichung als ,,mathematisches Modell‘ des untersuchten Vorganges bezeichnet 
werden. Beispiele für solche ,,mathematischen Modelle‘ sind die häufig ver- 
wandten Gleichungen = a—an?; »=a—fBn undn=a—-am—Pn. 
Sie sind im Gegensatz zu den oben erwähnten ,,reaktionskinetischen Modellen“ 
im allgemeinen nicht durch Elektronenübergänge zwischen verschiedenen 
Zuständen geschlossen deutbar, enthalten jedoch nur wenige Parameter und 
sind ohne Schwierigkeit geschlossen integrierbar. 


Solche „mathematischen Modelle‘‘ sind besonders bei der Diskussion von 
Leitfähigkeitsuntersuchungen herangezogen und häufig überbewertet worden. 
Die genauere Analyse gewisser experimenteller Ergebnisse lehrt, daß ihre 
einheitliche Beschreibung durch diese einfachen Modelle nicht möglich ist 
(vgl. Höhler?)). 

Bereits elementare Untersuchungen über das An- und Abklingen der 
elektrischen Leitfähigkeit bei Photoleitern nach unstetiger Anderung der 
optischen Anregung zeigen, daß bessere, „reaktionskinetische Modelle“ zur 
Diskussion des experimentellen Befundes herangezogen werden müssen. 

In dieser Arbeit soll als nichtstationärer Vorgang das Anklingen der Leit- 
fähigkeit nach dem Einschalten einer optischen Anregung genauer untersucht 
werden. Es soll zunächst gezeigt werden, welche Folgerungen für dieses 
Anklingen aus ganz allgemeinen reaktionskinetischen Modellen gezogen werden 
können. 


Dann soll das langzeitige Anklingen der Leitfähigkeit bei sehr schwacher 
optischer Anregung von CdS-Kristallen an einem speziellen reaktionskineti- 
schen Modell, bestehend aus einer Donatortermgruppe‘), mindestens einer 
Akzeptortermgruppe®) und dem Leitungsband genauer diskutiert werden 
(s. Abb. 2). 

In diesem Modell werden Elektronen durch die optische Anregung aus den 
Donatoren in das Leitungsband befördert. Von hier aus können sie mit den 
Übergangswahrscheinlichkeiten 8 bzw. y an Akzeptoren anhaften bzw. mit 
ionisierten Donatoren rekombinieren. Eine nicht verschwindende thermische 


8) G. Höhler, Ann. Physik 12, 379 (1953); 14, 426 (1954). 

*) Hier und im folgenden soll an Stelle einer optischen Anregung aus dem Valenzband 
nur eine solche aus Donatortermen angenommen werden (Ausläuferanregung), um von 
einer Berücksichtigung der Defektleitung absehen zu können. 

_ 5) Aus Gründen einer einheitlichen Terminologie benutzen wir im folgenden stets 
die Begriffe ,, Akzeptor“ und „Donator‘ ; dabei entsprechen die Donatoren den Aktivatoren, 
die Akzeptoren den Hafttermen der in der Lumineszenzphysik üblichen Bezeichnungweise. 
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Anregung gestattet zudem einen Ubergang von Elektronen aus Akzeptor- 


termen in das Leitungsband mit der Übergangswahrscheinlichkeit «. Diff 
Unter Berücksichtigung der Quasineutralitätsbedingung erhält man für 
dieses Modell folgendes Differentialgleichungssystem 
a+ah—Bn(H —h)—yn(n + h) Dal 
h=ßn(H—h)—ah, sinc 
sam 
wobei H die Konzentration der Akzeptorterme und Ah die Konzentration der geb 
in ihnen gefangenen Elektronen ist. a ist die Zahl der durch die optische An- kon 
regung pro sec und cm? ins Leitungsband beförderten Elektronen. | 
m, . 
und 
2. Allgemeine Überlegungen zum Anfangsteil von Anklingkurven von 
Belichtet man Photoleiter mit äußerst schwachen Lichtintensitäten, so oe 
stellt man häufig fest, daß sich die Leitfähigkeit auch längere Zeit nach dem rn 
Beginn der optischen Anregung noch nicht merklich ändert (dies ist insbe- 2 
A sondere dann der Fall, wenn y 
sich die Photoleiter zuvor lange Sat 
Zeit in Dunkelheit befanden : 
i oder ausgeheizt wurden). Erst di 
Einschalten der ganz allmahlich beginnt dann si 
Lichtanregung die Leitfähigkeit zu steigen, leit 
| wachst immer rascher mit der a 
| Zeit an und nähert sich nach 
j Durchlaufen eines Wende- = 
punktes einem stationären ma 
_t Wert (vgl. Abb. 1)*). Dieser eel 
to experimentelle Befund kann mit 
Abb. 1. Gemessenes Anklingen der Photoleitung den auf S. 11 genannten rein = 
bei sehr schwacher optischer Anregung. „mathematischen Modellen“ Zus 
(Schematisch) nicht beschrieben werden, weil : 
deren Lésungskurven bei physikalisch sinnvollen Anfangsbedingungen stets 2 
wendepunktsfrei sind. " 
Der in Abb. 1 deutlich erkennbare sehr allmähliche Leitfähigkeits- 
anstieg am Anfang des Anklingens legt die Vermutung nahe, daß im Punkt — 
t = t, die Steigung der Anklingkurve verschwindet, d.h. daß n(t,) = 0 gilt’). 
Wir werden im folgenden zeigen, daß diese einfache Annahme nicht berech- 
tigt ist. 
Nimmt man nämlich für Photoleiter an, daß die optische Anregung Ladungs- | 
träger aus nichtleitenden Zuständen direkt in den leitenden Zustand befördert 
(z. B. in das Leitungsband), so gilt für ganz allgemeine Modelle das folgende 
8) Vgl auch I. Broser u. R. Warminsky, Z. Physik 138, 340 (1952). _ 
?) Hier und im folgenden wird stets eine konzentrationsunabhängige Beweglichkeit, 
also Proportionalität zwischen der Konzentration der Stromträger und der Leitfähigkeit Ab 


angenommen (o = e nb). 
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Differentialgleichungssystem : 
nh = a—f (n, my, M,... m,) 

m, = f, (nm, m,...m,) (2) 

m, = f, (m, m,, . - m,)- 
Dabei ist n die Ladungsträgerkonzentration im leitenden Zustand; die m, 
sind Ladungsträgerkonzentrationen in nichtleitenden Zuständen. a ist die Ge- 
samtzahl der pro sec und em? durch die Anregung in den leitenden Zustand 
gebrachten Ladungsträger ; in den übrigen Gleichungen sind die entsprechenden 
konstanten Glieder in den Funktionen f, (n, m,,. . ., m,) einbegriffen. 

Alle Konzentrationen sind stetig in der Zeit, also auch die Funktion f(r, 

...m,), die den Ladungsträgeraustausch zwischen dem leitenden Zustand 
und allen nichtleitenden Zuständen beschreibt. Andert sich die Anregung z. B. 
vom Wert a= 0 bei t< t, unstetig auf den Wert a = a, beit >t), so muß 
auch » um den entsprechenden Betrag springen. War nun für t< t, die 
Elektronenkonzentration n stationär, d.h., war » = 0, so muß diese Ab- 
leitung bei t = t, auf einen endlichen Wert springen, und zwar in unserem 
Beispiel genau auf den Wert ay. 

Damit ergibt sich der für nichtstationäre Leitungsvorgänge wesentliche 

Satz: 
Erfolgt eine Anregung von Ladungsträgern aus nichtleitenden 
direktindenleitenden Zustand, so hat bei unstetiger Anderung der 
Anregung die Ableitung der Ladungsträgerkonzentration im 
leitenden Zustand hier ebenfalls eine Unstetigkeitsstelle. 

Eine weitere Aussage über die Anfangssteilheit der Anklingkurve gewinnt 
man leicht in einem etwas spezielleren Modell. In diesem wird analog dem 
obigen Modell eine direkte optische Anregung aus nichtleitenden Zuständen 
in den leitenden Zustand angenommen. Außerdem wird vorausgesetzt, daß 
ein direkter Ladungsträgeraustausch zwischen nichtleitenden Zuständen 
ausgeschlossen ist, d.h. daß ein solcher Austausch nur über den leitenden 
Zustand (z. B. das Leitungsband) erfolgt. A 

Ein Beispiel eines solchen Modells ist in Abb. 3 wiedergegeben ; es enthält 
eine tiefliegende Donatortermgruppe, beliebig viele Akzeptortermgruppen 
(in Verallgemeinerung des Modells von Abb. 2) sowie das Leitungsband. Die 


Leitungsband Leitungsband 


B a ay 4 a, 
h H, 


be 


Valenzband 


Abb. 2. Einfaches reaktionskinetisches Abb. 3. Reaktionskinetisches Modell mit 
Modell mehreren Akzeptortermgruppen 
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optische Anregung erfolgt aus den Donatoren in das Leitungsband (Ausläufer- 
anregung). Die Akzeptorterme können nur über das Leitungsband gefüllt 
und entleert werden. Das entsprechende Differentialgleichungssystem lautet: 


= — hy) — (3) 


hy = Ben (H, 


Es läßt sich nun leicht ne daß alle A, monoton mit der Zeit 
sind, und zwar im oben behandelten Fall des Anklingens monoton wachsend, 
d.h., es gilt stets A, > 0. 


Eine einfache Umformung der ersten Differentialgleichung des Systems (3) 
ergibt 


Da n, h,, h; und y niemals negativ sind, gilt ganz allgemein, daß n stets 
kleiner oder höchstens gleich dem Anregungsglied a ist: Es gilt für alle t: 
n(t) <a. 

Damit ergibt sich unter Benutzung des ersten Satzes ein zweiter für nicht- 
stationäre Leitungsvorgänge: 

Erfolgt eine Anregung von Ladungsträgern aus nichtleitenden 
direktindenleitenden Zustand undiisteindirekter Ladungsträger- 
austausch zwischen nichtleitenden Zuständen ausgeschlossen, 
sohatbeiunstetigerAnderungderAnregungdiezeitliche Ableitung 
der Stromträgerkonzentration im Zeitpunkt der Anregungs- 
änderung ihren maximalen Wert. 

Im Rahmen solcher Modellvorstellungen, die eine direkte Anregung von 
Ladungsträgern in den leitenden Zustand enthalten®), kann man also den in 
r Abb. 1 angegebenen Kurvenverlauf nur 
verstehen, indem man einen weiteren 


Wendepunkt bei kleinen Zeiten annimmt. 
n 


Die Steigung der Kurve n(t) kann dann 
| bei ¢ = t, tatsächlich ihren maximalen 
Wert a, haben, nimmt später jedoch 
bedeutend ab und beginnt erst nach 
längeren Zeiten zu wachsen, so daß die 
Kurve schließlich nach Durchlaufen eines 
nochmaligen Wendepunktes ihrem statio- 
nären Wert zustrebt: Im Bereich mitt- 


fo F 8) Vgl. Diplomarbeit H. Vogel, II. Phys. 
Abb. 4. Anklingen der Elektronen- Inst. der Humboldt-Universität, Berlin 1955. 
konzentration nach dem Einschalten 9) Im Gegensatz hierzu wäre zum Beispiel 
der optischen Anregung auch eine primäre Anregung in ein Excitonen- 

band denkbar. Die Excitonen tragen nicht 

zur Leitfähigkeit bei, können jedoch zum Beispiel nach thermischer Dissoziation 
freie Leitungselektronen liefern. Unter bestimmten Voraussetzungen wäre dann eine 
verschwindende Anfangssteigung von n(t) im Punkte t = t, möglich. 
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lerer Zeiten hat die Kurve ein ‚Plateau‘ (vgl. Abb. 4). Liegt dieses Plateau 
unterhalb der MeBgrenze, so ist es nicht nachweisbar und die MeBkurve hat 
das in Abb. 1 wiedergegebene Aussehen. 


Um die Vorstellung einer direkten Anregung von Elektronen in den lei- 
tenden Zustand bei Photoleitern vom Typ des CdS beibehalten zu können, 
muß man also zunächst nachweisen, daß es reaktionskinetische Modelle gibt, 
die solche ‚„‚Plateaukurven“ als Lösungskurven der entsprechenden Differential- 
gleichungssysteme enthalten. Außerdem sollten Bemühungen angestellt 
werden, durch Verfeinerung der Meßmethoden derartige ,,Plateaus“ direkt 
nachzuweisen, 


Beide Nachweise sind gelungen, worüber im folgenden berichtet wird. 
Es ist jedoch festzustellen, daß hiermit zwar ein notwendiges, aber keineswegs 
ein hinreichendes Kriterium für eine direkte Anregung von Elektronen in den 
leitenden Zustand gegeben ist. 


3. Diskussion eines reaktionskinetischen Modells mit einer Akzeptortermgruppe 


Es soll zunächst gezeigt werden, daß bereits das einfache Modell nach 
Abb. 2, für welches das Gleichungssystem (1) gilt, Lösungskurven enthält, die 
dem oben geforderten Verhalten entsprechen (Auftreten eines „Plateaus‘‘). 


Diese Lösungskurven werden in ihrer Gestalt wesentlich durch die Größe 
der im System (1) enthaltenen Parameter (der Übergangswahrscheinlichkeiten 
und Termkonzentrationen) bestimmt. 


Es wird sich nämlich zeigen, daß die Lösungskurven Wendepunkte gemäß 
Abb. 4 nur bei Erfüllung gewisser Relationen zeigen. 

Zur Vereinfachung der folgenden Betrachtungen wird zunächst ange- 
nommen, daß bei Beginn der Anklinguntersuchungen in den Akzeptortermen 
keine Elektronen mehr gefangen sind (das ist durch Ausheizen oder entspre- 
chend langes Warten im Dunkeln zu erreichen). 


Dann gilt für kleine Zeiten nach dem Einschalten der optischen Anregung 
stets, daß für die reaktionskinetischen Betrachtungen neben dem Anregungs- 
vorgang nur das Anhaften von Elektronen an den Akzeptoren wesentlich ist. 

Gegen diese Prozesse können also zunächst vernachlässigt werden: 

a) Die Reemission von Elektronen aus den Akzeptoren wegen des anfänglich geringen 
Füllungsgrades; diese Vernachlässigung gilt um so länger, je kleiner die Reemissionswahr- 
scheinlichkeit « ist. 

b) Die Rekombination von Elektronen mit ionisierten Donatoren, da die Rekom- 
bination quadratisch in ist; auch diese Vernachlässigung gilt um so länger, je kleiner die 
Rekombinationswahrscheinlichkeit y ist. 

Für „kleine Zeiten“ gilt dann an Stelle des Systems (1) das bedeutend 
einfachere System 

h=a—ßn(H —h), 
(la) 
h=ßn(H —h). 

Diese ,,kleinen Zeiten‘‘ werden dadurch begrenzt, daß die in (1a) vernach- 
lässigten, mit n und ¢ monoton wachsenden Glieder mit den übrigen vergleich- 
bar werden. 
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Ist der Bereich ‚kleiner Zeiten‘‘ groß genug, d.h. sind die Übergangs- 
wahrscheinlichkeiten x und y klein genug, so kann die Lösungskurve n({f) 
einen quasistationären Teil (ein „Plateau‘‘) enthalten!®): 

Die Elektronenkonzentration wächst zunächst mit der Steigung a; das 
Einsetzen des Haftvorganges begrenzt jedoch schnell (nach einer Zeit 
t, = 1/8 H) diesen Anstieg. Es kann sich also ein quasistationärer Bereich 
ausbilden (in (la) wird für % ~ 0 eine quasistationäre Elektronenkonzen- 
tration n ~ a/ß H erreicht, soweit h gegen H vernachlässigt wird). 

Bei Fortdauer der Einstrahlung werden die Akzeptoren immer mehr mit 
Elektronen gefüllt, und das für den Anhaftvorgang wesentliche Glied 
ön(H — h) nimmt monoton ab. — Würde keine Reemission aus Akzeptor- 
termen und keine Rekombination mit ionisierten Donatoren erfolgen, so 
wären alle Akzeptoren etwa nach einer Zeit tz = H/a gefüllt. — Es erfolgt 
dann ein weiteres lineares Anwachsen von n, das wie im Anfangsteil nahezu 
die Steigung a hat. Zugleich wachsen infolge der monotonen Zunahme der 
Haftelektronenkonzentration kh das Reemissionsglied x h und das Rekom- 
binationsglied yn (n + h). Es werden daher im folgenden zwei Fälle zu unter- 
scheiden sein, je nachdem, ob der Reemissions- oder der Rekombinationsvor- 
gang zuerst vergleichbar mit dem Anhaftvorgang wird"). 

Fall1: Das Reemissionsglied bestimmt den nun beginnenden Bereich 
„mittlerer Zeiten‘‘. 

Dann gilt 


n=a—ßn(H—h)+ah, 
kh=ßn(H—-h)—ah. 
Die Berücksichtigung der Reemission ergibt, daß der quasistationäre Teil der 
Anklingkurve schon endet, bevor die Akzeptorterme weitgehend gefüllt 


sind. Schon vor der Zeit r, wird die Elektronenkonzentration n wieder stärker 
ansteigen, nämlich zu einer Zeit 


(1b) 


H 
ap’ (4) 


die sich erwartungsgemäß um so weniger von r, unterscheidet, je kleiner « ist. 

Der Bereich ,,mittlerer Zeiten‘ wird dadurch nach oben begrenzt, daß das 
Rekombinationsglied y n(n + h) wesentlich wird. Die Elektronenkonzen- 
tration n wird jetzt nur noch weniger rasch anwachsen und sich schlieBlich 
dem wirklich stationären Endwert nähern. Dieser stationäre Endwert no 


berechnet sich aus dem System (1) für n = h = 0 und ergibt sich somit als 
Lésung der kubischen Gleichung 


Byrn? + (ByH (5) 


10) Diese Bedingungen sind also nur notwendig, nicht aber hinreichend für das Auf- 
treten eines Plateaus; im foıgenden werden auch die hinreichenden Bedingungen angegeben 
(s. Gl. (9) und (10)). 

11) Zur Entscheidung der Frage, ob das Reemissionsglied bei der vorliegenden Unter- 
suchungstemperatur wesentlich in Erscheinung tritt, kann zweckmäßigerweise die Ab- 
hängigkeit der stationären Elektronenkonzentration von der Anregungsintensität unter- 
sucht werden. Nur in den Bereichen, in denen das Reemissionsglied keine Rolle spielt, 
kann diese Abhängigkeit durch ein lineares bzw. ein Wurzelgesetz beschrieben werden 
(s. Gl. (6)). Im anderen Falle können auch von 0,5 und 1 bedeutend abweichende Ex- 
ponenten gefunden werden (vgl. Gl. (5)). 
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Fall 2: Das Rekombinationsglied bestimmt den Bereich ,,mittlerer Zeiten“. 
Dann gilt 
—h)—yn(n-+ Äh), 


h = Bn(H —h). ye 


Das Plateau kann sich in diesem Fall praktisch bis zur Zeit r, erstrecken. 

Bei weiterer Füllung der Akzeptorterme kann das Rekombinationsglied 
jedoch ein weniger rasches Anwachsen von n bewirken und schließlich zur 
Annäherung an einen stationären Endwert führen. Das Rekombinationsglied 
kann sogar, falls y groß genug ist!2), den Wiederanstieg der Elektronenkonzen- 
tration im Bereich ,,mittlerer Zeiten‘‘, (d. h. „nach Beendigung des Plateaus‘‘) 
völlig unterdrücken. Hier würde also kein stufenförmig abgesetztes Plateau in 
Erscheinung treten. 


In diesem Fall2 ergibt sich die stationäre Elektronenkonzentration als 
Lösung einer quadratischen Gleichung, die aus (1c) für n = h = 0 folgt: 


neo = (6) 


Unter der Voraussetzung, daß x und y klein genug sind, können also Lö- 
sungskurven n(t) auftreten, die Wendepunkte und sogar Plateaucharakter 
besitzen. Dabei wächst n(t) zunächst mit der Steigerung a an. Nach der Zeit 
t,~ 1/8 H kann n(t) zu einem quasistationären Plateau abbiegen. Nach 


einer Zeit 7, = = _ er steigt dann n(t) wieder stärker an und nähert sich 


schließlich seinem aus (5) bzw. (6) berechenbaren stationären Wert. 

Die oben angegebenen Voraussetzungen über x und y sind zwar notwendig, 
nicht jedoch hinreichend für das Auftreten eines Plateaus in n(t). 

Charakterisiert ist ein Plateau durch folgende Merkmale: 

1. t, ist bedeutend kleiner als r,; 

2. Die Steigung der Lösungskurve an ihrem ersten Wendepunkt ist be- 
deutend kleiner als die Anfangssteigung ; 

3. die Elektronenkonzentration am zweiten Wendepunkt ist bedeutend 
größer als am ersten Wendepunkt von n(t). 

Daraus ergibt sich eine weitere notwendige Bedingung über die Anre- 
gungsdichte a. Zusammen mit den vorgenannten Voraussetzungen eines 
kleinen « und y lassen sich diese Bedingungen, — wie durch eine einfache 
Rechnung gezeigt werden kann, — erfüllen, wenn 


n= zm<l, (8) 
(9) 


(10) 


2) Eine entsprechende Bedingung ist in Gl. (10) quantitativ formuliert. 
Ann. Physik. 6. Fo!ge, Bd. 17 


n (6) 
das 

Zeit 
reich 
IZEeN- 

mit 
lied 
Dtor- 
, 
folgt 
hezu 
<om- 
nter- 
reich 
| der 
füllt 
‚rker 

(4) 
ist. 
das 
zen- 
Blich 
Noe 
t als 

(5) 
Auf- 
eben 
nter- 
. Ab- und 
nter- 
rden 

.» 

g 


+18 Annalen der Physik. 6. Folge. Band 17. 1955 


gilt. Die Bedingungen (8), (9) und (10) sind notwendig und hinreichend für das 
Auftreten eines Plateaus in der Lösungskurve n(t) des Systems (1). 
Physikalisch sind diese Bedingungen prinzipiell für Kristalle, bei denen die 
Übergangswahrscheinlichkeit ß in die Akzeptoren größer ist als die Übergangs- 
wahrscheinlichkeit y in den Grundzustand, stets zu befriedigen. Bedingung (8) 
ist dann durch Wahl einer entsprechend geringen Intensität der optischen 
Anregung und Bedingung (9) 
N (kleines x) durch die Wahl 
i einer geniigend tiefen Unter- 
suchungstemperatur zu er- 
a füllen. 
| Umgekehrt läßt sich aus 
u einem Auftreten experimen- 
teller Anklingkurven mit 
Wendepunkten schließen, 
daß die Bedingungen (8), (9) 
und (10) erfüllt sind, d. h. 
daß insbesondere die Wahr- 
scheinlichkeit für den An- 
haftvorgang bedeutend 
größer als jene für die Re- 
kombination mit ionisierten 
Donatoren sein muß (Be- 
dingung (10)). 
Damit ist ein experimen- 
telles Kriterium zur Ent- 
Abb. 5. Theoretische Schar von Anklingkurven der scheidung der häufig inter- 


Elektronenkonzentration bei optischer Anregung. : a hd 
Scha licht essierenden Frage nach dem 
Verhältnis ß/y gegeben (z. B. 


darf beim Auftreten von Wendepunkten das Retrapping nicht vernach- 
lässigt werden). 

Abb. 5 zeigt nun eine schematische Schar von Lösungskurven des Systems 
(1) für den Fall 6 Sy. Scharparameter ist die Lichtintensität. 

Für sehr kleine Lichtintensitäten kann bei gegebenen Ubergangswahr- 
scheinlichkeiten und gegebener Akzeptorkonzentration H die Bedingung (9) 


oder die Bedingung (10), d.h. a >5 oder a > oH nicht mehr befriedigt 


werden. Hier nähert sich n(t) sehr langsam ohne Wendepunkt seinem statio- 
nären Wert. 


Für sehr große Lichtintensitäten ist die Bedingung (8) verletzt. Hier 
nähert sich n(t) sehr schnell, — ebenfalls ohne Wendepunkt, — seinem statio- 
nären Wert. 


Nur für mittlere Lichtintensitäten bilden sich typische Plateaukurven aus. 
Der Anfang des Plateaus 


(11) 


(siehe S. 17) ist praktisch unabhängig von der anregenden Lichtintensität. 
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Die Länge des Plateaus 


1 
nimmt umgekehrt proportional mit der Lichtintensität ab. Die Höhe des 


Plateaus 
© (12) 


(siehe S. 17) nimmt dagegen proportional der Lichtintensität zu. Das Ende 
aller Plateaus der Schar in Abb. 5 liegt also näherungsweise auf der Hyperbel 


Te (13) 


Für die Ausbildung eines besonders deutlichen Plateaus ist mithin eine . 
große Akzeptorkonzentration H sowie eine große Wahrscheinlichkeit ß für 
den Anhaftvorgang wesentlich. 


Umgekehrt lassen sich aus einer Meßkurve, die ein solches Plateau zeigt, 
durch Bestimmung der Plateaulänge 7, und der Plateauhöhe n, , recht gut 
die Werte für 6 und H entnehmen. i 


So ergibt sich fiir tiefe Temperaturen (« sehr klein; vgl. (4) und (5)) 


Hwa: (14) 


und 
1 


(15) 


p= 


Unter der Voraussetzung, daß nur eine Akzeptortermgruppe den Verlauf 
des Anklingens wesentlich bestimmt und insbesondere Einflußauf die Schaffung 
eines Plateaus in der Anklingkurve nimmt, ist hiermit eine Methode gegeben, 
ß und H für diese Termgruppe zu bestimmen. 


Weitergehende Untersuchungen an Modellen mit mehreren Akzeptorterm- 
gruppen unterschiedlicher energetischer Lage zeigen, daß die grundsätzlichen 
Überlegungen dieser Arbeit erhalten bleiben. Unter gewissen, durchaus reali- 
sierbaren Voraussetzungen über Konzentrationen und Übergangswahrschein- 
lichkeiten bleiben sogar die quantitativen Ergebnisse erhalten. Es treten 
dann unter Umständen mehrere verschieden ,,hohe‘‘ Plateaus in der Lösungs- 
kurve n(t) auf. Eine ausführliche Veröffentlichung über diese Untersuchungen 
ist in Vorbereitung. 


4. Experimentelle Ergebnisse an CdS-Einkristallen 


Die von uns untersuchten Kristalle waren durch Sublimation in einer 
H,—H,S-Atmosphäre aus einer hoch gereinigten, gefällten CdS-Ausgangs- 
substanz hergestellt. Sie waren nach einem von Muscheid und Buttler 
publizierten Verfahren randschichtarm mit Goldelektroden verspiegelt. 
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Die Kristalle standen in einer Vakuumapparatur mit einem kühlbaren 


Kupferblock in innigem Wärmekontakt. 


Die Messung erfolgte bei einem 


Vakuum von 10° Torr. Die Spannung wurde durch Goldkontakte zugeführt, 


welche auf die Verspiegelung 
federnd drückten (Meßfeldstärke 


etwa 50 V/cm). Der Strom wurde 
mit Hilfe eines Gleichstrom- 
verstärkers gemessen, der bei 


einer Zeitauflésung von etwa 
1 sec eine maximale Empfindlich- 
keit von etwa 2 - 10-12 A/Skt hatte, 

Die optische Anregung erfolgte 


mit einer Metallfadenlampe nach 
Monochromatisierung mit einem 
Schott-Filter (A = 547 mu) im 


L 


Ausläufer. Die Intensität wurde 
durch Graufilter variiert. Das 
Ein- und Ausschalten der Licht- 


10 107 
[sec] 


Abb. 6. Anklingkurven der Photoleitung bei 


sehr geringen Anregungsintensitäten 


10% 
t 


anregung wurde mit Hilfe eines 
Photoverschlusses durchgeführt 
(Verschlußkonstante < 1072 see). 

Die gesamte Meßapparatur 
wurde in einem lichtdichten Ge- 
häuse untergebracht, da bereits 


sehr kleine Lichtintensitäten — z. B. das Streulicht der Instrumenten- 
beleuchtung — ausreichen, um durch Akkumulierurg der wenigen hierdurch an- 


10% 
ter 
[sec] 
Abb. 7. Anklingkurve der Photoleitung 


bei sehr geringer Lichtintensität (MeB- 
temperatur: —100° C) 


geregten Elektronen in den Akzep- 
toren das Auftreten von Plateau- 
kurven zu verhindern. 

Vor der Aufnahme der Meßkurven 
wurden die Kristalle im Dunkeln 
etwa eine halbe Stunde lang bei etwa 
150° C ausgeheizt. Nach Abkühlung 
auf etwa —100°C ergaben sich die 
in Abb. 6 wiedergegebenen Kurven. 
Bei recht geringen‘ Anregungsinten- 
sitäten zeigt sich bei dem unter- 
suchten Kristall ein ausgeprägtes 
Plateau (Kurve 1). Wurde der Kristall 
jedoch nicht lange genug ausgeheizt, 
d.h. blieben die für die Ausbildung 
des Plateaus wesentlichen Akzeptor- 
terme bei Beginn der Messung noch 
teilweise gefüllt, so ergab sich 
Kurve 2, bei der dieses Plateau 
erwartungsgemäß weniger deutlich 
ausgebildet ist. Bei höheren Licht- 
intensitäten tritt kein Plateau mehr 
in Erscheinung (Kurve 3). 
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Das Ende des Plateaus der Kurve 1 liegt bei etwa 100 sec. Die Elek- 
tronenkonzentration, die diesem Plateau entspricht, betrug etwa 10! cm-*. 
Bei einer Anregungsintensität, die etwa 1013 Elementarakten pro cm? und sec 
entspricht, ergibt sich für den untersuchten Kristall eine Akzeptorkonzen- 
tration von H = 1015 cm? und eine Übergangswahrscheinlichkeit von 
B » 10713 cm~* sec-!, also durchaus plausible Werte. 


Bei entsprechender Wahl der Anregungsbedingungen (Temperatur- und 
Liehtintensität) ist es durchaus möglich, auch mehr als ein Plateau zu finden, 
was auf eine entsprechende Wirksamkeit mehrerer energetisch verschieden 
liegender Akzeptortermgruppen hinweist. In Abb. 7 ist eine solche Meßkurve 
für den gleichen Kristall gezeichnet. Es sind deutlich zwei Plateaus bei etwa 
30 und bei 200 sec zu erkennen. 


Für eine ausführliche Diskussion dieser Kurve mit mehreren Plateaus 
muß auf die in Vorbereitung befindliche zweite Arbeit verwiesen werden. 


Bei den hier angegebenen Messungen handelt es sich zunächst um orien- 
tierende Untersuchungen, die lediglich zur Bestimmung der Größenordnung 
des zuvor beschriebenen Effektes sowie zu einer mehr qualitativen Bestätigung 
der theoretischen Überlegungen angestellt wurden. 


Genauere experimentelle Untersuchungen, insbesondere mit weitgehend 
monochromatischer Anregung"), sind in Vorbereitung. 


5. Zusammenfassung 


Es wurde gezeigt, daß bei direkter Anregung von Elektronen in leitende 
Zustände (Band) die Anklingkurve der Leitfähigkeit ganz allgemein — unab- 
hängig vom weiteren Reaktionsmechanismus — mit endlicher Steigung be- 
ginnen müssen. Diese Anfangssteigung ist die maximale überhaupt auftretende 
Steigung und gleich der Zahl der Anregungsprozesse pro cm? und sec. 


Es wurde weiter gezeigt, daß das Auftreten der experimentell bereits 
bekannten Wendepunkte im Anklingen der Photoleitung bei geringer An- 
regungsintensität wesentlich durch Akzeptoren bedingt ist und daß man aus 
der Lage dieser Wendepunkte Aussagen über die Akzeptorkonzentration 
sowie die Übergangswahrscheinlichkeit für den Haftprozeß machen kann. 


Die Parameterbereiche (insbesondere Lichtintensität und Temperatur), 
für die solehe Anklingkurven mit Wendepunkten bzw. ausgeprägten Plateaus 
auftreten, wurden diskutiert. 


Abschließend wurde eine Reihe orientierender Messungen an CdS-Ein- 
kristallen beschrieben, die Anklingkurven mit solchen Wendepunkten bzw. 
mit ausgeprägten Plateaus liefern. Die Größenordnungen der aus diesen Meß- 


13) Die zwischen Ausläufer- und Grundgitteranregung um viele Größenordnungen 
veränderliche Absorptionskonstante geht entscheidend in die Auswertung der Ergebnisse 
ein (Bestimmung des Anregungsgliedes a); sie erfordert bei Ausläuferanregung eine sehr 
weitgehende Unterdrückung kurzwelligen Streulichtes. 
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kurven berechneten Werte von H und ß stimmen mit den aus anderen Uber. 
legungen (z. B. Glowkurven) zu erwartenden Werten überein. 


Berlin, Laboratorium für die Physik des elektrischen Durchschlages im 
Institut für Strahlungsquellen der Deutschen Akademie der Wissenschaften 
und II. Physikalisches Institut der Humboldt-Universität. 


Bei der Redaktion eingegangen am 7. April 1955. 
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Die Amplitude des primären Echos 
bei der Hahnschen Spin-Echo-Methode 


Von Harry Pfeifer 
Mit 3 Abbildungen 


Inhaltsiibersicht 


Es wird die Amplitude des primären Echos als Funktion der Resonanz- 
abweichung berechnet. Die entstehenden Kurven werden insbesondere im 
Hinblick auf die Messung der transversalen Relaxationszeit diskutiert und 
mit experimentellen Ergebnissen verglichen. Abschließend wird der Einfluß 
der Feldinhomogenität auf die Echoamplitude näherungsweise untersucht. 


Die Hahnsche Spin-Echo-Methode!) eignet sich besonders zur Messung 
der transversalen Relaxationszeit 7,?). Legt man nämlich an die zu unter- 
suchende Probe Materie, die sich in einem magnetischen Gleichfeld H, befindet, 
nacheinander zwei gleiche hochfrequente Rechteckimpulse der Frequenz 
y H,/2 x (y = gyromagnetisches Verhältnis der betreffenden Atomkerne), so 
entsteht nach dem zweiten Impuls eine resultierende Magnetisierung, das 
sogenannte Echo, dessen Amplitude als Funktion des Abstands der beiden 
Impulse in einfacher Weise 7, ergibt. Bei der Aufnahme einer derartigen MeB- 
reihe muß der zeitliche Abstand zwischen zwei aufeinanderfolgenden Impuls- 
paaren wenigstens größer als die thermische Relaxationszeit 7, sein*). Diese 
liegt bei den meisten Substanzen in der Größenordnung einiger Sekunden, 
so daß die Resonanzbedingung w = y H, über einen längeren Zeitraum erfüllt 
sein muß. Da aber sowohl das Magnetfeld H, als auch die eingestrahlte Frequenz 
zeitlichen Schwankungen unterworfen sind, ist es von größtem Interesse, 
die Abhängigkeit der Echoamplitude von der Resonanzabweichung zu kennen, 
um hieraus entstehende Fehler bei der Messung der transversalen Relaxa- 
tionszeit so klein wie möglich zu halten. Die hierfür entwickelten Gleichungen 
gestatten außerdem, den Einfluß der Feldinhomogenität auf die Echoamplitude 
näherungsweise zu untersuchen. 


Setzen wir voraus, daß die Impulsdauer t, klein ist gegen T,, T, und den 
Abstand r der beiden Impulse, und wählen wir den Zeitnullpunkt zu Beginn 
des ersten Impulses, so folgt aus den Blochschen Gleichungen?) für die 


1) E. L. Hahn, Physic. Rev. 80, 580 (1950). 

*) N. Bloembergen, E. M. Purcell, R. V. Pound, Physic. Rev. 78, 679 (1948). 
8) H. Pfeifer, U. Dolega, H. Winkler, Ann. Physik 15, 246 (1955). 

‘) F. Bloch, Physic. Rev. 70, 460 (1946). 
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Magnetisierung in Richtung der magnetischen Komponente der Hochfrequenz 
impulse: 

M, = mW cos wt — vy sin wt, (1) 
wobei 


Ble \ (vy Terme mit 
Ue=2) = 2H, me (sin 2 ( B ) B cos (nt (y H, — w) + 9,) 
?2r 
Btu\/yHı\® | Terme mit 
= H, (sin B ty) (sin 2 B ) + cos (nt (y H, — 0) + Pn) 


ß®= (y H, — + y2H} n = 1,2. (2) 


H, ist dabei die halbe Amplitude der magnetischen Komponente der Hoch» 
frequenzimpulse und x, die statische Kernsuszeptibilität. In allen praktischen 
Fällen wird das Magnetfeld H, über die Probe nicht konstant sein und damit 
die Magnetisierung sich als Mittel der für die verschiedenen Felder resul 
tierenden Magnetisierungen ergeben. Legen wir eine Gau ßsche Feldverteilung 
zugrunde, so wird (für 73< r) 


-d(y H, —y H,) (3) 


3t _ #2 


-d(y H,—y H,). 


V2ln2/T% ist dabei die halbe Halbwertbreite der Feldverteilungsfunktion, 
deren Zentrum bei H, liegt. Setzen wir voraus, daß y H, 73 >1,t,/T3 ist, 
so liefert die Integration 


\3 * 


y A, 


(4) 


ar 
= Ay % 
1 


(sin)? Hj + @*?t,) 


wobei w* = y H, — w die Resonanzabweichung darstellt. Die Echoamplitude 
selbst ergibt sich zu 


= Hy % (sin? = Hi + o*:) a 


1 


+ w* 
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Wie bekannt, hat das Echo im Resonanzfalle Maxima fiir y A, t,, = 2 2/3 +2kn 
und 47/3 +2k2;k=0,1,2... Das sind also die Fälle, die bei einer 
Messung der transversalen Relaxationszeit praktischerweise in Frage kommen. 
In Abb. 1 ist die Echoamplitude als Funktion der Resonanzabweichung für 
die ersten vier Maxima dargestellt. 


Echoamplitude 
(willkürliche Einheiten ) 


Y = 7 


YA, 
Abb.1. Die Amplitude des primären Echos als Funktion der Resonanzabweichung für 
yH,t„=2n/3, 42/38, 20/3 + 2a, + 22 


Man sieht, daß für k = 0 nur ein Echomaximum nämlich im Resonanzfalle 
auftritt, während für k = 1 Nebenmaxima auftreten. Für k=2,3...ent- 
stehen weitere Nebenmaxima, wie man leicht nachprüfen kann. Um für eine 


a) 


Abb. 2. Das primäre Echo als Funktion der Resonanzabweichung. a) Für y H, t, = 42/3; 
b) Fir y H,t, = 22/3 + 22. (Zur Verdeutlichung wurde in Abb. 2b die Echoein- 
hiillende nachgezeichnet) 


gegebene Abweichung von der Resonanz die Anderung der Echoamplitude 
möglichst klein zu halten, muß man also mit möglichst kräftigen Hochfre- 
quenzimpulsen arbeiten (y H, groß) und die Impulsdauer t,, so wählen, daß 
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y H,t, =“* wird. Die Abb. 2 zeigt für die beiden Fälle y Hy t, = 42/8 
und y H,t, = 22/3 + 2a das primäre Echo als Funktion der Resonanz- 
abweichung. Durch Variation und Messung der Impulsbreite t,, kann man auf 
diese Weise aus der Kurvenform bequem die Impulsamplitude 2 H, bestimmen, 
was meßtechnisch wesentlich einfacher ist als eine direkte Messung von 2 H,. 
Abweichungen von der Rechteckform der Hochfrequenzimpulse und Ver- 
stimmungen der Verstärkerkreise machen sich besonders bei größeren Impuls- 
breiten bemerkbar. Die Nebenmaxima besitzen dann unterschiedliche Ampli- 
tude, fehlen mitunter ganz oder es treten noch zusätzliche Nebenmaxima auf. 

Die Gln. (3) gestatten, die Echoamplitude im Resonanzfalle näherungs- 
weise als Funktion der Inhomogenität (73) zu berechnen, und zwar ergibt sich 


50% 
fehlen 
kleine 


u = 0, 


da uw unsymmetrisch in y H,—w ist und für » unter der Voraussetzung 
Hy 


1 1 1 ty 
yi+ (=) sin (y Hy ty + are 
1 
sin (2 y H,t,, are tg (= = (6) 
Dieses Ergebnis geht für t,/y H, Tf? = 0 in das von Hahn?) angegebene 
Resultat über und ge- 
stattet unter den obigen 
Voraussetzungen, den 
Einflußder Feldinhomo- 
genität auf die Echo- 
amplitude zu beurteilen. 
al In Abb. 3 ist die Echo- 
amplitude als Funktion 
von t,/H,y Tz? dar- 
13 gestellt. 

Man sieht, die maxi- 
male Echoamplitude 
wird für t,/yH, Ty? ~ 

Y 0,2 erhalten, jedoch 
verläuft im betrach- 
10 - > teten Bereich das 
tw Maximum sehr flach. 

yy T°? Zum Beispiel macht 

Abb.3. Die Echoamplitude im Resonanzfalle als Funktion eine Änderung der Feld- 
von t,,/y H, T¥? verteilungsbreite um 


5) E. L. Hahn, l.c. 


Echoamplitude 
(willkürliche Einheiten ) 


= 


H. Pfeifer: Amplitude des primären Echos bei der Hahnschen Spin-Echo-Methode 27 


Amplitudenänderung aus. Es ist also lediglich emp- 


50% weniger als 5% 
nzurichten, daßt,/yH, 


fehlenswert, die experimentellen Bedingungen so ei 
kleiner als 1 ist. 


Leipzig, Physikalisches Institut der Karl-Marx-Universität. 


Bei der Redaktion eingegangen am 15. April 1955. 
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Ein Gerät zur automatischen Messung dieses 
des Wolterschen Farbabstandes') Licht 
Beleu 
funkt 
Mit 8 Abbildungen funkt 
RB. 
Inhaltsübersicht 
Es wird ein kleines Spektral-Schnellphotometer mit automatischem daher 
Rechengerät beschrieben, das den Abstand zweier Farben (oder Beleuch- auf e 
tungen) im Sinne der Wolterschen Definition?) anzeigt. Das Gerät dient @ Licht 
insbesondere zur Prüfung künstlicher Lichtquellen (z. B. Leuchtstofflampen) sprec. 
hinsichtlich ihrer spektralen Abweichungen von einer „Sollbeleuchtung“ Ergel 
(z. B. Tageslicht). Der angezeigte Abstand ist ein Maß für die Ersetzbarkeit zweie 
zweier Beleuchtungen bezüglich ihrer unterschiedlichen Wirkung auf den fasseı 
Farbeindruck der beleuchteten Körper. D 


Farb 
vorlie 
dense 
punk 
liche 


Von Hans-Walter Bodmann 


Einleitung 


Die begriffliche Fassung des Wortes „Farbe“ muß in der Physik das Ziel 
haben, eine objektive Eigenschaft einer Strahlung, eines Körpers usw. zu AR 
kennzeichnen. Das wird im Hinblick auf die Bedürfnisse der physikalischen griff? 
Meßtechnik gefordert, die Farben als Hilfsmittel zur Messung physikalischer stanc 
Größen heranzieht. Als Beispiel sei das zweidimensionale Farbschlieren- strah 
verfahren?) erwähnt, das die Lichtablenkungen nach Richtung und Betrag Diffe 
durch Farbton und Sättigung erkennen läßt. Nach Wolter‘) verstehen wir werd 
daher unter dem Wort Farbe die spektrale Energieverteilung einer Licht- ziert 
strahlung oder das spektrale Remissionsvermögen eines Körpers (damit ist wird 
allgemein ein „„Farboperator‘‘ gemeint, der jeder auffallenden Strahlung eine setzu 
neue, von dem Körper ausgehende Strahlung zuordnet). D 


In der Praxis ergibt sich häufig die Aufgabe, die Unterschiede zweier Verg 
Farben geeignet zu charakterisieren. Zur quantitativen Fassung dieser Frage siert 
benötigt man einen Begriff für den „Abstand‘‘ zweier Farben. Im folgenden Farb 
betrachten wir speziell den Farbabstand von Lichtstrahlungen, die zur Be- in eit 


1) Dissertation Kiel 1954; vorgetragen auf der Tagung des Verbandes der Deutschen 
Physikalischen Gesellschaften in Hamburg am 18. September 1954. 


2) H.W. Bodmann, Ann. Physik (6) 12, 348 (1953). verte 

3) H. Wolter, Ann. Physik (6) 8, 1 (1950). 

*) H. Wolter, „Ziele und Wege der Farboperatorenlehre‘, Vortrag vor dem Deut- brau 
schen Fachnormenausschuß ‚Farbe‘ in Göttingen am 26. September 1951. defir 
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leuchtung von Körpern dienen. Damit soll die Frage nach der Ersetzbarkeit 
von Lichtquellen erfaßt werden. 


Als Ausgangspunkt betrachten wir einen Körper, der mit der Strahlung 
einer Lichtquelle beleuchtet werde. Das von dem Körper reflektierte Licht 
falle in unser Auge und verursache einen Farbeindruck (z. B. „rot‘‘). An Hand 
dieses Vorganges bezeichnen wir die spektrale Intensitätsverteilung J, der 
Lichtquelle im folgenden als Beleuchtungsfunktion oder auch einfach als 
Beleuchtung, das spektrale Reflexionsvermögen R, des Körpers als Remissions- 
funktion und schließlich die spektrale Augenempfindlichkeit V, als Bewertungs- 
funktion. Der von uns registrierte Farbeindruck wird dann durch das Produkt 
I;- R,: V, bestimmt. Ersetzen wir nun die Lichtquelle durch eine andere 
mit der Beleuchtungsfunktion /f, so ändert sich im allgemeinen auch der 
Farbeindruck, den ein Beschauer von dem Körper erhält. Wir betrachten 
daher zwei Lichtquellen als identisch, wenn ihre Funktionen J, und I} bis 
auf einen konstanten Faktor übereinstimmen. Inwieweit z. B. technische 
Lichtquellen das Tageslicht ersetzen, kann im einzelnen nur aus den ent- 
sprechenden Beleuchtungsfunktionen entnommen werden. Es liegt nahe, das 
Ergebnis dieses Vergleiches durch eine geeignete Maßzahl für die Ersetzbarkeit 
zweier Lichtquellen auszudrücken. Dieser ‚Abstand‘ muß dann möglichst um- 
fassend die Unterschiede der Beleuchtungsfunktionen kennzeichnen. 


Der Vergleich der entsprechenden Farbtemperaturen oder der Punkte im 
Farbendreieck ist gerade bei multiplikativer Farbenmischung, wie sie hier 
vorliegt, unzulänglich. So werden z. B. sehr verschiedene Lichtstrahlungen durch 
denselben Punkt im Farbendreieck dargestellt. Daher können die Farb- 
punkte zweier Lichtquellen den Abstand null haben, obwohl sie unterschied- 
liche Farbeindrücke eines Körpers liefern. 


H. Wolter‘) hat im allgemeineren Zusammenhang einen Abstandsbe- 
griff für „„Farben‘‘ angegeben, dem die Hilbertsche Definition für den Ab- 
stand zweier Funktionen zugrunde liegt. Im wesentlichen werden zwei Licht- 
strahlungen in der Weise miteinander verglichen, daß man über die quadrierte 
Differenz der spektralen Energien integriert und geeignet normiert. Dabei . 
werden die Verteilungsfunktionen mit einer Bewertungsfunktion multipli- 
ziert (bei Beleuchtungen z. B. der spektralen Augenempfindlichkeit). Hier 
wird insbesondere der Abstand nur dann null, wenn die spektrale Zusammen- 
setzung der Strahlungen übereinstimmt. 


Die Bedeutung dieses Abstandsbegriffes erstreckt sich allgemein auf den 
Vergleich zweier Lichtstrahlungen als ‚Farben‘. Insbesondere charakteri- 
siert er die Ersetzbarkeit zweier Beleuchtungen, soweit Veränderungen des 
Farbeindruckes der beleuchteten Körper in Betracht kommen. Das zeigte sich 
in einer Arbeit), in der aus spektroskopischen Messungen der Abstand einiger 
technischer Lichtquellen vom Tageslicht bestimmt wurde. Für die meßtech- 
nische Anwendung sind derartige numerische Rechnungen sehr zeitraubend, 
abgesehen von der vorher notwendigen Messung der spektralen Energiever- 
verteilungen. Es soll daher die Aufgabe dieser Arbeit sein, die Bestimmung 
des Wolterschen Farbabstandes zu automatisieren. Zu diesem Zweck 
braucht man ein Gerät, das für vorgegebene Lichtquellen das in der Abstands- 
definition enthaltene Rechenprogramm selbständig ausführt. 
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Absehnitt 1 


Den Farbabstand zweier Lichtstrahlungen mit den relativen spektralen 
Intensitätsverteilungen J, (A) und J, (A) schreiben wir zunächst in der Form 


{@, q I, (A) — G, 1, (A)}? da 
@ = 5 Mi (1) 


wenn G, eine spektrale Gewichtsfunktion bedeutet, die den Eigenschaften des 
maßgebenden Empfängers Rechnung trägt°). Bei Anwendung dieses Ab- 
standsbegriffes auf Beleuchtungen z. B. verstehen wir unter G, eine Be 
wertungsfunktion des Auges®). Der variable Faktor q ist mit Rücksicht auf 
die optimale Anpa-sung der relativen Verteilungsfunktionen angebracht. 
Das Minimum wird bestimmt durch 


13 (A) dA 


(A) da 


Damit wird 


oo 
1, (A) - I, (A) dA 


@=1 ——,0>d>1. (3) 


Varwa-fanaa 
v 0 


Den Abstand d = 0 haben dabei nur proportionale Intensitätsverteilungen 
I (A). Der größtmögliche Abstand d = 1 besteht z. B. zwischen solchen Licht- 
strahlungen, die keinen Spektralbereich gemeinsam haben (allgemein zwischen 
orthogonalen Funktionen G, J;). 


Es zeigt sich nun, daß diese Abstandsdefinition für die praktische An- 
wendung noch erweitert werden muß. Betrachten wir z. B. als Beleuchtungen 
I, (A) und I, (A) zwei benachbarte, sich nicht überdeckende Spektrallinien, so 
wird ihr Abstand nach (3) gleich eins. Andert sich aber die Reflexion aller 
beleuchteten Körper in diesem Spektralbereich sehr wenig, so wird ihr Farb- 
eindruck praktisch nicht geändert, wenn die Linien genügend nahe beiein- 
ander liegen. Wenn (3) eine sinnvolle Aussage über die Ersetzbarkeit zweier 
Beleuchtungen ergeben soll, müssen wir bei der Messung der Beleuchtungs- 
spektren eine genügend große Spaltweite wählen, die von der Flankensteil- 
heit in den Remissionsfunktionen der beleuchteten Körper bestimmt wird. 
Allgemein sind alle Funktionen 7, im Spektrum noch zu ‚verschmieren“. 
Dadurch läßt sich der Abstandsbegriff an die ,,Kérperfarben“ anpassen. 


5) Wir betrachten hier nur den Farbabstand zweier Lichtstrahlungen in bezug auf 
einen Empfänger, d.h. G,(A) = @, (A). 

6) Bei Messungen mit Farbfilmen z. B. wäre eine durch deren Eigenschaften bedingte 
Gewichtsfunktion zu wählen. 


standi 
Grund 
Vekto: 


findlich 


30 

Die n 

oo 

Sie ist 

| Vi @ (A) da- (A) da in Fra 

terscl 

4 

oo Die D 

m. (2) Zu: 

J gemeir 

De 

die für 

kann 

Die 

begriff 

Farbeı 

Du 

einer s 

Die 

Darste 

indem 

Darin 


H.-W. Bodmann: Ein Gerät zur automatischen Messung des Wolterschen Farbabstandes 31 


Die mathematische Formulierung führt zu einer Faltung der bewerteten 
Verteilungen G, : J, mit einer „Spaltweitenfunktion‘‘ ®; (4’ — A): 


f(A) = —A) (4) 


Darin ist ®, eine Rechteckfunktion, Gaußsche Fehlerfunktion oder dgl. 
Sie ist so zu definieren, daß ihre Breite den Flankensteilheiten der überhaupt 
in Frage kommenden Körperfarben angemessen ist’). Somit wird der Wol- 
tersche Farbabstand endgültig 


oo (oo oo 
0 Ww 0 


oo (wo 2 00 (wo 2 
V f ®, (a —a)G, 2’) aa da f f I, da 
0 
(5) 


Die Definition (3) ist hierin als Spezialfall enthalten. Sie geht aus (5) hervor, 
wenn man als ®, die Diracsche 6-Funktion wählt. 

Zusammenfassend erwähnen wir noch einmal die Eigenschaften der all- 
gemeinen Definition (5): 

Der Farbabstand enthält zwei noch wählbare Parameter (Funktionen), 
die für die jeweilige Fragestellung von der Erfahrung bestimmt werden. (Diese 
kann psychologisch-physiologischer und physikalischer Natur sein.) 

Die Definition (5) besitzt alle Eigenschaften des Hilbertschen Abstands- 
begriffes für Funktionen. Insbesondere wird d= 0 nur für proportionale 
Farben J;. 

Durch die Normierung wird (5) von den absoluten Intensitäten sowie von 
einer speziellen Normierung der Funktionen G, und ®, unabhängig. 

Schließlich ist die Wellenlängenskala keineswegs ausgezeichnet. 


Die geometrische Deutung dieses Abstandsbegriffes schließt sich an die 
Darstellbarkeit von Funktionen durch Ortsvektoren im Hilbert-Raum an, 
indem wir für (4) schreiben 


f(A) = Px (A), A, 


Darin sind die fy abzählbar viele Koordinaten nach der orthogonalen, voll- 
ständigen Basis der px(A) (z. B. das System der Fourier-Funktionen). Auf 
Grund der Eigenschaften der px(A) können wir den Quotienten in (5) dann in 
Vektorschreibweise formulieren : 


y 
L 


?) Bei Beleuchtungsfunktionen kann durch ®, auch die spektral verschiedene Emp- 
findlichkeit des Auges fiir Farbunterschiede beriicksichtigt werden. 
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Das entspricht dem ,,Kosinus des Winkels‘ zwischen /,(A) und f,(A). Geome 
trisch ist also der Abstand (5) allein durch den Winkelabstand, der die Funk 
tionen /, (A) und f, (A) repräsentierenden Ortsvektoren im Hilbert-Raum 
bestimmt. 

Wir spezialisieren nun die allgemeine Form der Definition (5) auf den 
Farbabstand zweier Beleuchtungen. Zunächst machen wir eine Festsetzung 
über @,. Um eine definierte Funktion vor Augen zu haben, wählen wir für 
G, probeweise die internationale spektrale Hellempfindlichkeit V, des Auges, 
Damit ergibt sich ein für das Auge maßgeblicher Integrationsbereich sowie eine 
angemessene Bewertung der spektralen Unterschiede zweier Beleuchtungen, 
Weiterhin haben wir aus den Absorptionsgebieten der Körperfarben eine ge- 
eignete Spaltweitenfunktion ®; festzulegen entsprechend den in Frage. kom- 
menden Flankensteilheiten. Für die vorliegende Arbeit haben wir bei der 
Messung der Beleuchtungsspektren versuchsweise eine Spaltbildweite von 
150 A im mittleren Teil des sichtbaren Spektralbereichs gewählt. Das ent 
spricht etwa den in diesem Teil des Spektrums auftretenden Absorptions 
flanken der praktisch in Frage kommenden Körperfarben. 

Nach dieser Wahl von @, und ®; schreiben wir die Abstandsdefinition in 
der Form, an die sich das im folgenden beschriebene Gerät anschließt. Nach (2) 
gilt zunächst für das Minimum von (1) 


Gi Gain. 12 (A) dA = 63 12 (A) da. 


Mit der Bezeichnung qyın. * I} (A) = I, (A) und J, (A) = I(A) wird die Defini- 
tion (1) dargestellt durch 


Ip (2) —G, ar 


2 


O<d<l: (7) 


oo 
(ayaa 
u 


wenn wir als Normierungsbedingung fordern 


(Darin sei z. B. /, (A) die Beleuchtungsfunktion des Tageslichtes, / (A) die einer 
Glühlampe.) 

An Hand von (7) und (8) ergibt sich das Programm für eine automatische 
Abstandsbestimmung. Die spektrale Differenz der zu vergleichenden Licht 
strahlungen läßt sich mit Hilfe eines Zweistrahl-Spektralapparates dureh 
photoelektrische Empfänger in Differenzschaltung registrieren. Diese Diffe 
renz-Photospannungen sind dann zu quadrieren und schließlich über alk 
sichtbaren Spektralbereiche durch ein Element genügend großer Zeitkor 
stante zu mitteln. Dieser Mittelwert stimmt bis auf einen unwesentlichen 
Faktor, der bei der Normierung herausfällt, mit dem Zähler von (7) übereit, 
wenn als Apparatefunktion die spektrale Augenempfindlichkeit realisiert ist 
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Die Normierungsbedingung (8) läßt sich durch Abschwächen einer der Strah- 
lungen immer erfüllen. Durch abwechselndes Abdecken jeweils eines Strahlen- 
ganges kann man die Gleichheit der betreffenden Integrale prüfen und den 
Nenner von (7) ablesen. Die Ausdehnung des Empfängers in Richtung des 
Spektrums bestimmt dabei das gleichzeitig zur Messung kommende mittlere 
Spektralintervall, d. h. die mittlere effektive Spaltbildweite. 


Abschnitt 2 


In der Abb. 1 ist der prinzipielle Aufbau des Gerätes in halbschematischer 
Form dargestellt. Die zu vergleichenden Lichtstrahlungen fallen jeweils auf 
eine Hälfte des Spaltes eines Spektralapparates. Unmittelbar hinter dem 
Eintrittsspalt (in Abb. 1 vor dem Spalt gezeichnet) befindet sich ein Synchron- 
Motor mit einer Sektorscheibe, die das durch den Spalt eintretende Licht zum 
Zwecke der nachfolgenden Meßverstärkung mit 1500 Hz moduliert. Da die 
erforderliche Auflösung sehr gering ist, entsprechend einer mittleren Spalt- 
bildweite von 150 A, genügt als Spektralzerleger ein Glasprisma mit kleiner 
Dispersion. Die Basislänge beträgt hier 3,5 cm bei einem Öffnungsverhältnis 
des Spektralapparates von 1:10. 
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DIFFERENTIAL-\ 
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VERSTÄRKER $37) ANZEIGE 


Quadrierung und 
Integration 


spektr Differenz Normierung 
Abb. 1. Die wesentlichen Elemente des Gerätes und ihre Funktionen 


Das spektral zerlegte Licht beider Strahlungen wird durch einen Kipp- 
spiegel über einen photoelektrischen Differenzempfänger geführt, der sich am 
Bildort des Spektrums befindet. Dieser Empfänger besteht aus zwei gegen- 
einandergeschalteten Selen-Photoelementen, die von einem mittleren Spektral- 
bereich von 150 Ä ausgeleuchtet werden. Ihre spektrale Empfindlichkeit 
gestattet eine weitgehende Angleichung an die Augenempfindlichkeit. Hier 
erfolgt die Anpassung durch geeignete Wahl der spektralen Abtastgeschwin- 
digkeiten. Die Winkelgeschwindigkeiten des Kippspiegels werden erzeugt 
durch einen gleichförmig umlaufenden Kurvenkörper. Sein Rand bewegt 
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Abb. 2. Differentialphotoelement. Die 
Breite betrigt 1mm. Das bedeutet am 
Bildort des Spektrums einen mittleren 
Wellenlängenbereich von 150A. Die 
Kriimmung entspricht dem Bild des 
Spaltes. In der Mitte ist die lichtemp- 
findliche Schicht unterbrochen, so daß 
zwei Elemente mit gemeinsamer Gegen- 
elektrode entstehen. Diese stimmen in 
ihrer spektralen Empfindlichkeit sehr 
genau überein 
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einen mit der Spiegelachse fest verbundenen Schwenkarm. Die Kurvenscheibe 
ist nun so geformt, daß hinsichtlich der Integration die gesamte Apparate 
funktion gleich der spektralen Augenempfindlichkeit wird. Das sichtbar 
Spektralgebiet wird 10mal in der Sekunde abgetastet. 


Uber einem linearen Potentiometer von 5 k Q wird die Differenz der Phote- 
spannungen abgenommen und einem Meßverstärker zugeführt. Der Schleifer 
des Potentiometers liegt an der gemeinsamen Gegenelektrode beider Elemente, 
Diese erhalten durch Verändern der Schleiferstellung verschiedene Arbeits- 
widerstände, die es gestatten, die absoluten Intensitäten beider Lichtquellen 
im Sinne der Normierungsbedingung aufeinander abzustimmen. Der Ausgang 
des Verstärkers ist auf ein Thermokreuz angepaßt, das die verstärkte Differenz- 
spannung quadratisch gleichrichtet. Die Trägheit des Thermokreuzes dient 
gleichzeitig zur Integration über alle Spektralintervalle. Mit Hilfe eine 


. Form des verwendeten Empfängers dar- 


strahlungsstärke von 2 - 10-5 W/cm?. Das bedeutet maximal eine Beleuch- 
tungsstärke des Empfängers von 10? Lux (Abb. 3). Bei dem verwendeten 
Spektralapparat liefern technische Lichtquellen im allgemeinen geringere Be- 
strahlungsstärken. Intensivere Strahlungen können durch Graugläser oder 
geschwärzte Siebe immer entsprechend geschwächt werden. Nach Abb. 3 


empfindlichen Gleichspannungsinstru- 
mentes kann man schließlich Zähler 
und Nenner von (7) direkt ablesen. 

In den nächsten Abschnitten folgen 
Einzelheiten über den Empfänger, über 
die Meßverstärkung sowie über die An- 
gleichung an die Augenempfindlichkeit. 


Abschnitt 3 


Der photoelektrische Empfänger hat 
hier die Aufgabe, die spektrale Differenz 
zweier Lichtstrahlungen zu messen. 
Dabei soll seine spektrale Empfindlich- 
keit möglichst gut mitder Augenempfind- 
lichkeit übereinstimmen. Diese Bedin- 
gung wird am besten erfüllt von Selen- 
Photoelementen. In der Abb. 2 ist die 


gestellt, der aus einem kommerziellen 
Selen-Photoelement hergestellt wurde. 

Der Photostrom (und damit der 
Spannungsabfall über dem Arbeitswider- 
stand) eines Selen-Photoelementes ist für 
nicht zu große Bestrahlungsstärken pro- 
portional der auffallenden Lichtinten- 
sität, sofern der Arbeitswiderstand ge- 
nügend klein ist. Hier verläuft die 
Kennlinie für einen Arbeitswiderstand 
von 5kQ noch ganz im Proportionali- 
tätsbereich bei einer maximalen Be- 
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steht jedenfalls ein relativer Intensitätsmeßbereich von zwei Zehnerpotenzen 
zur Verfügung. 

Eine absolute Empfindlichkeitsmessung im Bereich um 5460 Ä ergab 
0,1 Ampere pro Watt oder etwa 2 - 10-* Ampere pro Lux und cm? wirksamer 
Fläche. Die Nachweisgrenze liegt hier bei 5 - 10? Watt (Signal-Rauschver- 
hältnis 10:1) oder mit einer wirksamen Fläche von 0,06 cm? bei etwa 0,5 Lux. 


Intensitöt Intensität Intensität 
— — 


0 


0% 50% 700% 700% 50% 100% 


Abb. 3. Linearität von Photospannung und Lichtintensität für drei Spektralbereiche, 

gemessen bei maximalen Beleuchtungsstärken bis 100 Lux. Ordinate ist die Spannung am 

Verstärkerausgang. Zur Intensitätsschwächung dienten geschwärzte Siebe mit photo- 

metrisch gemessenen Schwächungsfaktoren. Die Abweichungen liegen innerhalb der 

Meßgenauigkeit dieser Faktoren. Die Rauschspannung betrug 0,15 V. Der Maximal- 
ausschlag wurde jeweils auf 15 V eingeregelt. 


Spektrale Empfindlichkeitsunterschiede zwischen beiden Elementen ließen 
sich nicht nachweisen. Zur Prüfung wurden beide Elemente von einem Spek- 
trum gleichmäßig ausgeleuchtet. Dabei ergab sich bei durchlaufenden Spek- 
trum und voller Verstärkung außer der reinen Rauschspannung kein Signal. 


Abschnitt 4 


Die von dem Differentialphotoelement gelieferten Spannungen sind hier 
von der Größenordnung 10% und 10-3 Volt. Daher schließt sich an den 
Empfänger ein mehrstufiger Resonanzverstärker an (Abb. 4). Die Gesamt- 
verstärkung beträgt 10°. 

Als Bandbreite wurde 300 Hz gewählt derart, daß Modulationsschwan- 
kungen vom Netz her praktisch keinen Verstärkungsabfall bedingen. Hin- 
sichtlich der Abtastgeschwindigkeiten des Spektrums ist hier die Zeitkonstante 
des Verstärkers belanglos, da über lange Zeiten gemittelt wird. Die Linearität 
der Verstärkung ist in Abb. 3 enthalten. Das Eigenrauschen des Verstärkers 
beträgt einige 10-7 Volt in guter Übereinstimmung mit dem berechneten 
Wert 2 - 10-7 Volt. 
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Abb. 4. Photospannungsverstärker für Wechsellicht von 1500 Hz 


H.-W. Bodmann: Ein Gerät zur automatischen Messung des Wolterschen Farbabstandes 37 


Am Ausgang des Verstärkers erfolgt die Quadrierung der Photospannung 
mit Hilfe eines Vakuum-Thermokreuzes (TH 5), dessen Heizfadenwiderstand 
von 1,1Q über einen Aus- 
gangstransformator an den 
Innenwiderstand der End- Br, 


röhre ECC40 angepaßt ist. 
Die EMK des Thermo- 


elementes ist im Aus- 
steuerungsbereich des Ver- 
stärkers (bis 15 Volt) sehr 
genau proportional dem 
Quadrat der Ausgangsspan- 
nung (Abb. 5). 

Die Zeitkonstante des 
Thermokreuzes beträgt etwa Ug? 
3sec. In dieser Zeit wird - 
das sichtbare Spektrum 1 


1 10 "100 Voir? 


ers ers. =>)» m 
Galvanometer (RG4) ver- 


wendet, das von dem niederohmigen Thermoelement sehr stark bedämpft 
wird. 


Abschnitt 5 


Zu der von Optik und Empfänger bedingten Apparatefunktion kommt hin- 
sichtlich der Integration als weitere Gewichtsverteilung die Abtastzeit pro 
Spektralintervall. Durch geeignete Wahl der spektralen Abtastgeschwindig- 
keiten ergibt sich somit eine kontinuierliche Angleichungsmöglichkeit an die 
spektrale Augenempfindlichkeit. Der Vorteil dieser Methode gegenüber Ver- 
wendung von Filtern liegt in der optimalen Ausnutzung der Empfängeremp- 
findlichkeit. 


In Abb. 6 ist der verwendete Abtastmechanismus dargestellt. Der Schwenk- 
arm des Kippspiegels wird mit Federkraft an eine Kurvenscheibe gedrückt, 
die von einem Synchronmotor in Umdrehung versetzt wird. Damit der Motor 
möglichst gleichmäßig belastet wird, befindet sich außer dem Spiegelarm 
ein leerlaufender Schwenkarm an der Kurvenscheibe.. Zwischén beiden ist 
eine (regelbare) Feder aufgehängt. Ihre Rückstellkraft ist groß gegen die Träg- 
heitskräfte der Schwenkeinrichtung bei der Maximalbeschleunigung. Zur 
Stabilisierung des Umlaufs läuft auf der Achse der Kurvenscheibe außerdem 
ein Schwungrad, an das eine Keilriemenuntersetzung vom Antriebsmotor her 
angreift. Das Schwungrad trägt eine Gradeinteilung, die bei gleichférmigem 
Umlauf mit gleichen Zeitabständen identisch ist (t-Skala). Die Gleichförmig- 
keit des Umlaufs wurde stroboskopisch geprüft. 


Zur Bestimmung der Form der Kurvenscheibe, die die richtigen spektralen 
Winkelgeschwindigkeiten des Kippspiegels liefert, mißt man zunächst die 
Abhängigkeit A = A (r), wenn r den variablen Abstand vom Drehpunkt zum 
Rand der Kurvenscheibe bezeichnet. Diese bleibt unabhängig von der Form 


> 
> 
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der Scheibe erhalten, so daß man aus der noch zu bestimmenden Funktion nisch tı 
A = Alt) auf die gesuchte Kurvenform r = r (t) schließen kann. Die zeitliche @ mit Hil 
Abtastung 4 = A(t) erhält man nun durch Vergleich der schrittweise ge- = 
ten hie 
silber- 
Empfanger Bogenli 
Zwische 
interpo 
trum / 
2 eines Sp 
Vo j ram-Gl 
Tempe! 
(2700 - 
Kurvenscheibe Mit de 
Kippspiege! Absorp 
ram si 
spektra 
Auges 
V (A) - 
beschri 
einen F 
in Abb 
In 
Abb. 6. Die Wobbeleinrichtung des Gerätes. (Die A- und t-Skala dient nur zur Eichung) das = 
gestellt 


der spe 
messenen Photospannungen U (A) eines bekannten Spektrums /° (A) mit der „n ea 


berechneten Verteilung V(A)-J° (A). Diese werden von vornherein nicht schen : 
proportional sein. Wir setzen daher eine z 
dt durch 

V (A)- 1° (A) = K-U (A)- (A), K= const. Abb. 8 


stellun 


Schwungrad Vom Synchr Motor 


Daraus folgt durch Integration 
Die 
dA = t (A), t,, K = const. (WE Geräte 
tremer 
durch 
Dieses Integral läßt sich numerisch auswerten, wobei A von A, bis A,, zwei # neten 
Grenzwellenlängen im violetten und roten Teil des sichtbaren Spektralbereichs, ®$ durch | 
läuft. Innerhalb dieses Bereiches soll die Angleichung an die Augenempfind- @ Werte 
lichkeit erfolgen. Außerhalb dieses Intervalles geben wir uns zwei Wellenlängen 8 stände 
A, und A, vor, für die t = 0° und t = 180° festgesetzt wird: großer 
Absta 
„<A P sung n 
0° < t,< t,< 180°. (10) 100-W 
ladun; 
Damit sind die Konstanten K und t, festgelegt. Die Gebiete links und 100-W 
rechts von A, bzw. A, werden für die Umkehr des Kippspiegels unter mecha- @ Stift 


U (A) 


@ 
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nisch tragbaren Beschleunigungen benötigt. Aus (9) und (10) läßt sich dann 
mit Hilfe von A = A (r) die Kurvenscheibe punktweise berechnen. 


Zur Wellenlängeneichung dien- 
ten hier die Linien einer Queck- 
silber-Entladungslampe und die 
Bogenlinien des Lithiums. Die 
Zwischenwerte wurden rechnerisch 
interpoliert. Als bekanntes Spek- 
trum /° (A) wurde die Strahlung 
eines spannungsstabilisierten Wolf- 
ram-Glühfadens benutzt, dessen 
Temperatur pyrometrisch zu 
(2700 + 30)° K gemessen wurde. 
Mit den tabellierten Werten des 
Absorptionsvermégens von Wolf- 
ram sowie der internationalen 
spektralen Hellempfindlichkeit des 
Auges läßt sich die Verteilung 
V (A) - I° (A) berechnen. Das oben 
beschriebene Verfahren lieferte 
einen Kurvenkörper, dessen Form 
in Abb. 7 dargestellt ist. 

In der folgenden Abb. 8 ist 
das Ergebnis der Angleichung dar- 
gestellt. Die genaue Anpassung 
der spektralen Gewichtsverteilung 
an die Augenempfindlichkeit zwi- 
schen 4,2 und 6,8 - 103 A erfordert 
eine zeitliche Abtastung A(t), die 
durch die gestrichelte Kurve in 
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Abb. 7. Der Kurvenkörper für die Angleichung 
der Apparatefunktion an die spektrale Augen- 
empfindlichkeit. Für einige t-Werte sind die 
zugehörigen Wellenlängen am Rande ver- 
merkt. Die zweite, zur ersten symmetrische 
Hälfte des Kurvenkörpers gehört zum Rück- 
lauf des Spektrums. Der Durchmesser der 
Symmetrielinie beträgt 60 mm. Die Kurven- 
scheibe erzeugt einen Winkelhub des Kipp- 
spiegels von 5°. Die radiale Toleranz beträgt 
0,01 mm 


Abb. 8 wiedergegeben ist. Die ausgezogene Kurve wurde nach der Her- 
stellung der Kurvenscheibe gemessen. 


Abschnitt 6 

Die Meßgenauigkeit des 
Gerätes wurde in zwei ex- 
tremen Fällen ermittelt 
durch Vergleich der berech- 
neten Abstände mit den 
durch das Gerät angezeigten 
Werten. Als Kontrollab- 
stände wurden ein sehr 
großer und ein sehr kleiner 
Abstand gewählt. Die Mes- 
sung mit dem Gerät ergab für 


100-W-Glühlampe/Hg - Ent- 
ladungslampe d=0,67 +0,03, 


7 
107A 


100-W-Glihlampe/Nernst- pp. 


t 


0° 45° 90° 135° 180° 
8. Die Abhängigkeit 2 = 4 (t). Gestrichelte 


Stift d= 0,075 + 0,005. Kurve berechnet, ausgezogene Kurve gemessen 
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Für die Berechnung dieser Abstände wurden die relativen spektralen Energie 
verteilungen von Glühlampe, Hg-Entladungslampe und N ernst-Stift mit 
dem Gerät als Spektralphotometer gemessen. Zu diesem Zweck wurde die 
Abtastzeit des Spektrums mit Hilfe einer großen Untersetzung soweit herauf- 
gesetzt, daß die spektralen Photospannungen direkt von einem Registrier 
galvanometer aufgezeichnet werden konnten. Als Normalstrahler diente 
dabei wieder der spannungsstabilisierte Wolfram-Glühfaden (7 = 2700° Kj. 
Mit den so bestimmten Energieverteilungen berechnen sich nach Gl. (3) von 
Seite 30 die Abstände 


100-W-Glühlampe/Hg-Entladungslampe d = 0,60, 
100-W-Glühlampe/Nernst-Stift . d = 0,083. 


Der Vergleich mit den obigen Werten zeigt Abweichungen von etwa 10%, 
Dieser Fehler enthält die Abweichung der Apparatefunktion von der Augen 
empfindlichkeit V,, die für die Rechnung benutzt wurde. Im Hinblick auf 
die nur vorläufige Festlegung der Bewertungsfunktion (vgl. S. 31) erhält 
also der Meßfehler eine geringere Bedeutung. 

Als erste Anwendung des Gerätes wurden die Abstände einiger moderner 
Lampentypen von Tageslicht gemessen. Diese Messungen wurden bei gleich- 
mäßig grau bedecktem Himmel gemacht. Die Ergebnisse unserer Messungen 
sind in der folgenden Tabelle zusammengestellt. 


Lampentype Abstand vom Tageslicht 


Gutam: 
Philips Tageslicht de Luxe (Leuchtstofflampe) . . . 
Philips Warmton de Luxe (Leuchtstofflampe) . . . 
Philips Weiß de Luxe (Leuchtstofflampe) . . . 
Philips Tageslicht (gefilt. Glühlampe) . . . 


OO Or 00 OO 


dad 


Das Variationsintervall für den Farbabstand zweier Lichtquellen er 
streckt sich von d = 0 (identisch) bis d = 1. Als Anhaltspunkt für die Be 
urteilung der obigen Werte sei noch angeführt, daß der Farbabstand von 
Sonnenlicht zu Tageslicht in der Größenordnung 0,1 liegt?). 

Es zeigt sich, daß die Abstände der untersuchten Typen nur wenig variieren, 
Im ganzen sind alle aufgeführten Lampen vergleichbar hinsichtlich ihrer Ver- 
wendbarkeit als Tageslichtersatz. Die hier aufgeführten Leuchtstofflampen 
zeigen wesentliche Verbesserungen gegenüber einem früher untersuchten 
Modell, das verglichen mit einer normalen Glühlampe einen sehr viel ungün- 
stigeren Abstand vom Tageslicht ergab. Bei diesem Modell älterer Bauart 
wurde d = 0,42 gemessen?). 

Für die Prüfung serienmäßiger Lampentypen braucht man eine konstante 
und definierte Normalbeleuchtung. Insbesondere wird ein Tageslichtnormal 
erforderlich, um unabhängig von Wetterbedingungen Messungen ausführen und 
miteinander vergleichen zu können. Im Zusammenhang damit sei auf die 
DIN-Beleuchtung®) zur Messung von Filmempfindlichkeiten hingewiesen. 


8) DIN-Normblatt Nr. 4512. 
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Als Lichtquelle dient dort eine genormte Metallfadenlampe, die in Verbindung 
mit einem speziellen Flissigkeitsfilter®) eine weitgehende Annäherung an das 
Spektrum von Sonnenlicht ergibt. 


Zusammenfassung 


Eine automatische Messung des Wolterschen Farbabstandes in der Form 
(7) und (8) (Seite 32) läßt sich durch Kombination eines Zweistrahl-Spektral- 
photometers mit einfachen elektrischen und mechanischen Elementen reali- 
sieren : 

Die Differenzbildung erfolgt mit Hilfe eines Differential-Photoelementes 
als Empfänger. Durch geeignete Wahl der spektralen Abtastgeschwindigkeiten 
werden beide Spektren mit der Augenempfindlichkeit bewertet. Die Differenz- 
Photospannungen werden verstärkt und in einem Thermokreuz quadratisch 
gleichgerichtet. Auf Grund seiner Trägheit integriert das Thermokreuz in 
Verbindung mit einem stark bedämpften Anzeige-Galvanometer gleichzeitig 
über alle Spektralintervalle. Zur Normierung nach Gl. (8) dienen regelbare 
Arbeitswiderstände der Photoelemente. 

Eine Kontrolle der Meßgenauigkeit des Gerätes ergab einen mittleren 
Fehler von etwa 10%. 

Durch Wahl der Spaltbildweite und durch Auswechseln einer Kurven- 
scheibe können die Parameter ®, und G, des Farbabstandes (vgl. Def. (5), 
S. 31) verändert werden. 


Dem Direktor des Instituts, Herrn Prof. Dr. Lochte-Holtgreven, 
dessen wertvolle Unterstützung die Durchführung dieser Arbeit ermöglichte, 
und Herrn Doz. Dr. Wolter, der die Arbeit anregte und durch stete Anteil- 
nahme förderte, möchte der Verfasser an dieser Stelle seinen aufrichtigen 
Dank aussprechen. 


®) R. Davis u. K. G. Gibson, Publ. Bureau of Stand Nr. 114. 


Kiel, Institut fir Experimentalphysik der Universitat. 


Bei der Redaktion eingegangen am 23. April 1955. 
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Zur Messung der Dielektrizitatskonstanten 
von Flüssigkeiten im Dezimeterwellenbereich 


Von Theobald Jäkel 
Mit 5 Abbildungen 


Inhaltsübersicht 


Es wird eine Methode zur Bestimmung der komplexen Dielektrizität 
konstanten von Flüssigkeiten im Dezimeterwellengebiet beschrieben, die eim 
Auswertung der Meßdaten ohne Verwendung von Hyperbelfunktionen ode 
von entsprechenden Näherungen gestattet. Die Messung erfolgt in einer 
konzentrischen Lecherleitung, die teilweise mit der Flüssigkeit gefüllt ist 
Bei konstanter Wellenlänge ergeben sich Realteil und Imaginärteil der Di 
elektrizitätskonstanten explizit aus der Anderung der Resonanzwerte der Meß 
leitung mit Variation der Probendicke. 

Nach einem Überblick über die wesentlichen Merkmale der bekannten 
Meßverfahren wird die Theorie der neuen Methode entwickelt. Dann folgt 
eine Beschreibung der Meßapparatur und der Messung, und schließlich werden 
einige Meßergebnisse für Alkohollösungen mitgeteilt. 


1. Einleitung 


Die Bestimmung der komplexen Dielektrizitätskonstanten & = e’ — je 
ist unerläßlich, wenn man die Theorien über die dielektrische Relaxation, über 
das innere Feld und die Assoziation in polaren Flüssigkeiten prüfen und 
klären will. Die sichersten Schlüsse können gezogen werden, wenn e’ und € 
im ganzen Bereich der anomalen Dispersion mit der Wellenlänge und de 
Temperatur als Parameter gemessen wurden. Die entsprechenden Reihen 
messungen im Höchstfrequenzgebiet sind langwierig, weil der experimentelk 
Aufwand mit zunehmender Frequenz wächst, aber auch weil die Berechnung 
der Dielektrizitätskonstanten aus den Meßwerten sehr umständlich ist. B& 
den komplizierten Auswerteverfahren ist eine laufende Kontrolle der Met 
sungen nicht möglich, außerdem läßt sich der Fehler des. Ergebnisses schlecht 
abschätzen. 

Wir beschränken uns auf das Dezimeterwellengebiet, also auf Verfahren, 
bei denen im allgemeinen Lecherleitungen verwendet werden. Die bekannte 
Meßmethoden schließen oft die Bestimmung der magnetischen Stoffkonstanten 
mit ein, was z. B. für die Untersuchung von Halbleitern wichtig ist. Da jedoeb 
die komplexe Permeabilität der polaren Flüssigkeiten fast immer gleich eim 
gesetzt werden kann, wollen wir die Methoden nur in bezug auf den rein di 
elektrischen Fall betrachten. 
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Allgemein wird folgender Weg beschritten: Es wird der Reflexionsfaktor 
einer Probe bestimmt, die in eine konzentrische Leitung eingeführt und mit 
einem bekannten Widerstand abgeschlossen ist. Man kann hierfür das Ab- 
tastverfahren benützen, das in der heute verwendeten Form von Knol und 
Strutt!) angegeben wurde, oder das Resonanzverfahren. Beim ersteren ergibt 
min 


sich der Betrag des Reflexionsfaktors f = o - e/¥ aus dem Verhältnis . 


max 
der Spannungsverteilung auf der Meßleitung, während man die Phase y 
aus der Entfernung des ersten Spannungsminimums vom Eingang der Probe 
berechnet. Eine wichtige Voraussetzung ist dabei, daß die Abtastsonde die 
Feldverteilung nicht merklich stört und daß ihre Eintauchtiefe konstant ist. 
Im Resonanzverfahren bestimmt man den Phasenwinkel y aus der Resonanz- 
länge der Meßleitung, den Betrag o aus der Halbwertbreite der Resonanz- 
kurve. Die wesentliche Voraussetzung ist hier die einwandfreie Funktion der 
verschiebbaren Kurzschlußbrücke. Schwan?) hat nachgewiesen, daß die 
Resonanzmethode dem Abtastverfahren gleichwertig ist und unter gewissen 
Bedingungen den Vorzug verdient. 


Aus dem Reflexionsfaktor kann der auf den Wellenwiderstand der Meß- 
leitung bezogene Eingangswiderstand der Probe errechnet werden, aus dem 
man schließlich die Dielektrizitätskonstante gewinnt, wenn die Probenlänge s 
und der Abschlußwiderstand bekannt sind. Meist arbeitet man mit offener oder 
kurzgeschlossener Probenleitung. Aus der Vierpoltheorie ergeben sich in diesen 
beiden Fällen transzendente, komplexe Gleichungen für &, die die Form 
tv = 4- tg 4 bzw. tw = 4- ctg4 haben. Darin ist tw eine komplexe Größe, die 
durch die Meßwerte gegeben ist, während in 4 neben der Wellenlänge A in Luft 
und der Probendicke s das gesuchte é enthalten ist. Die Lösung der Gleichungen 
erfordert Tabellen oder Diagramme der Funktionen; eine graphische Dar- 
stellung für weite Bereiche findet man bei Untermann?). 


izität 
lie eim 
n oder 
1 einer 
lt ist 
ler Di. 
r Meß. 


ınnten 


1 folgt 


verden 


Schwan hat in einer weiteren Arbeit?) die Auswerteverfahren ausführlich 
diskutiert und dabei jeweils den Umfang der notwendigen Rechenarbeit ange- 
geben. Er leitete unter der Voraussetzung kleiner Probendicken Näherungs- 
gleichungen für e’ und e’”’ ab, bei deren Verwendung man einen Fehler von 
höchstens 1% bzw. 3% macht. Es sind 22 Rechen- und Tabellenoperationen 
notwendig, wenn man die Auswertung mit Hilfe der Funktion 3-tg3 vornimmt, 
durch die Näherung erniedrigt sich diese Zahl auf 14. 

In manchen Fällen ist es angebracht, die Probenleitung mit einem belie- 
bigen, bekannten Blindwiderstand abzuschließen. Jungnickel und Fal- 
kenhagen) geben die Theorie zu einer Anordnung, bei der sich hinter der 
Probe eine zusätzliche kurzgeschlossene Leitung von variabler Länge be- 
findet. Diese Abschlußleitung ermöglicht es auch, die Fälle Kurzschluß und 


1) K.S.Knolu. J. O. Strutt, Physica IX, 577 (1942). 

*) H. Schwan, Ann. Physik (6) 5, 253 (1950). 

5) G. Untermann, Z. angew. Physik 2, 233 (1950). (Untermann bestimmt den 
Reflexionsfaktor nach dem Abtastverfahren. Er kommt zu einem falschen Ergebnis fiir 
die Phase, weil er zur Vereinfachung 2 hinzufügt, was einer Vergrößerung der Entfernung 
des ersten Spannungsminimums vom Probeneingang um 4/4 und nicht — wie angegeben — 
um 4/2 entspricht. Außerdem enthalten verschiedene Endformeln Rechenfehler.) 

*) H. Schwan, Ann. Physik (6) 5, 287 (1950). 
5) H. Jungnickel u. H. Falkenhagen, Ann. Physik (6) 9, 341 (1951). 
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Leerlauf einwandfrei zu verwirklichen, indem man ihre Länge gleich > 5 baw.4 Vor 
wählt, denn bei offenem Ende stören z. B. die Strahlung und die Streckapesiill PB 
Das Meßverfahren von Huber®) zeichnet sich dadurch aus, daß keine ver Es 
schiebbaren Kurzschlußkolben und Sonden verwendet werden, die in elek. 
trischer und mechanischer Hinsicht oft Fehlerquellen darstellen. Die Ab Z\ 
stimmung auf Resonanz oder die Abtastung erfolgt durch Verändern d ay 
Probendicke, der Betrag des Reflexionsfaktors wird aus der Maximalspan- 
nung und aus der Leerlaufspannung an der Sonde ermittelt, wobei letztere in = 
einer besonderen Messung bestimmt wird. Wenn nicht an einer kurzgeschloss- IL 
nen Probe in einer senkrecht stehenden Meßleitung gemessen wird, muß die 
Probe in eine Meßkammer aus Isoliermaterial gebracht werden, die die Le. sI 
tungshomogenität zusätzlich stört. Schwan hat in einer dritten Arbeit‘) G- 
untersucht, wie Halterungen am Ende von Lecherleitungen den Eingang 
widerstand der Probe beeinflussen und entsprechende Korrekturformeln ap. G- 
gegeben. Von Esau und Prösdorf wurde in einer kurzen Originalmitteilung‘) @= 
über ein Meßkammerverfahren berichtet, das nur sehr geringe Substanz k 
mengen erfordert und eine Auswertung mit Hilfe einfacher Näherungsglei Z 
chungen ermöglicht. Das Hauptproblem liegt hier in der geeigneten Formung 
der Kammer aus einem Material, das sich in der Meßflüssigkeit nicht löst. > ak 
Bei fast allen Methoden werden die Eigenverluste der Meßleitung und die _ 
Kopplungsdämpfungen vernachlässigt. Jungnickel und Falkenhagen‘) Wir be 
leiten eine Näherung für die Eigendämpfungskonstante ab, mit deren Hille die Vi 
experimentell die Korrektur bestimmt werden kann, die am Reflexions- 
faktor anzubringen ist. Gewöhnlich wird sie notwendig, wenn an Flüssig 
keiten mit geringen dielektrischen Verlusten im unteren Dezimeterwellen- 
gebiet gemessen wird oder wenn die Probendicke sehr klein ist. EN 
2. Theorie des neuen Verfahrens — 
Es ist naheliegend, ein Verfahren zu suchen, bei dem in der Auswertung 
keine transzendente Gleichung auftritt und außerdem der Umweg über den 
Reflexionsfaktor nicht notwendig ist. 
Wir setzen eine Meßanordnung voraus, die in Abb. 1 schematisch dar E5 Wi 
gestellt ist. Wir unterscheiden zwei Teile: Länge 
a) Die Meßleitung, die auf der einen Seite mit einem verschiebbar Differ 
Kolben K kurzgeschlossen und auf der anderen Seite mit dem Eingang der Term 
Probenleitung abgeschlossen ist. In der Kurzschlußbrücke befinden sich Kopp 
lungsschleifen für den Generator 
und den Indikator. 
=D b) Die Probenleitung, die di w;. 
é gleichen geometrischen Querdimen 
Sp sionen wie die Meßleitung besitzt 
K und die Meßflüssigkeit enthält. Sie 
Abb. 1. Schematische Darstellung der Meß- ist mit dem komplexen Leitwert 6, 
anordnung abgeschlossen. 
*) O. Huber, Z. angew. Physik 6, 9 (1954). Er 
?) H. Schwan, Ann. Physik (6) 5, 268 (1950). he, 
8) A. Esau u. D. Prösdorf, Naturwiss. 41, 208 (1954). richte 
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Vorausgesetzt sei fernerhin, daß die Ohmschen Verluste in der Leitung und 
die Kopplungsdämpfung vernachlässigbar sind. 


Es werden die folgenden Bezeichnungen verwendet: 


Z Wellenwiderstand der Meßleitung, 
A Meßwellenlänge in Luft, 
“= Phasenmaß der Meßleitung, 


| Länge der Meßleitung bei Resonanz, 

s Länge der Probenleitung, 

&=@ + j® Eingangsleitwert der Probenleitung, 
6-Z=g+ jp reduzierter Eingangsleitwert, 


e= e' — je” Dielektrizitätskonstante der Probe, 


Z=2. 1_ Wellenwiderstand der Probenleitung, 


ye 
y= jaVe Fortpflanzungskonstante der Probenleitung. 


Wir betrachten zunächst nur die Probenleitung, für deren Eingangsleitwert 
die Vierpoltheorie®) die Beziehung ergibt: 


26. + Tg (y 8) 


= (1) 


Beachten wir Z’ = und y = jx //éund lösen (1) nach ‘- auf, so erhalten 


ve 


wir: 
26, VéIg(jaVé s)—ZG 


VE 


Es wird nun angenommen, daß der AbschluBleitwert @, bei Variation der 
Länge s konstant bleibt und wir fragen nach der Anderung von Z& mit s. 
Differenzieren wir (2) nach s, so verschwinden im Zähler der rechten Seite die 
Terme mit den Hyperbelfunktionen und es folgt: 


(2) 


a .1 d(Z6) 
(3) 
Wir setzen Z& = g + jg, spalten (3) in Realteil und Imaginärteil auf und 
erhalten : 


(4) 


é (5) 


Fe i B. R. Feldtkeller, Einführung in die Vierpoltheorie der elektrischen Nach- 
tichtentechnik, Leipzig (S. Hirzel) 1948. 
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e’ und e’”’ können also aus dem Verlauf des reduzierten Leitwertes bei 
Variation von s berechnet werden, ohne daß der Abschlußleitwert ©, bekannt 
zu sein braucht. Die Länge der Probenleitung geht auch nicht in die Be 
ziehungen ein, was vorteilhaft ist, weil oft die Querschnittsebene nicht genau 
20 festgelegt werden kann, in de 
g / die ungestörte Feldverteilung 


t endet. 
es Für die Diskussion von (4) 


und (5) sind in Abb. 2 Realteil 
und Imaginärteil von Z & gegen 
s aufgetragen worden. Die an 
genommenen Zahlenwerte ent- 
sprechen einer Probenleitung, 
die mit einer Streukapazität be 
lastet ist. 

Der Verlauf von @ weist se 
fort auf eine Möglichkeit hin, die 
Gleichungen für e’ und e”’ zı 
vereinfachen, indem man be 
einer Probenlänge mißt, für die 
\ y, y = 0 wird, also bei s,, oder 4 

Der Realteil g besitzt in der 
Nähe der Nullstellen von p Ex 


10 


tremwerte. Aus (5) folgt mit 


Abb. 2. Reduzierter Eingangsleitwert Z@— — 0 Ohne weiteres, daß das 
9+j9 in Abhängigkeit von der Länge s bei Minimum von g bei einem s< sm 

2 0,33pF liegt. Aus meBtechnischen Grün- 

den ist im vorliegenden Falk 

die Messung an der Nullstelle vorzuziehen, der ein kleiner Wirkleitwert 
zugeordnet ist. 


Man könnte nun mit Hilfe einer vorgeschalteten MeBleitung nach der Ab 
tastmethode oder nach dem Resonanzverfahren den Verlauf von Z & in einem 
kleinen Bereich um s, aufnehmen und die Dielektrizitätskonstante aus 


erhalten. Dabei ist aber noch die Berechnung des reduzierten Eingangsleit- 
wertes notwendig. Es soll jetzt gezeigt werden, daß man e’ und e’’ ohne diesen 
Umweg aus den Meßwerten gewinnen kann. 


Djakov!P) gab eine Resonanzmethode zur Messung komplexer Wider- 
stände an, die die entsprechende Umformung gestattet. Die MeBleitung 
gleicht der in Abb. 1 dargestellten Anordnung, ihre Eigenverluste seien zu 
nächst vernachlässigt. Durch Umrechnung der Endformeln von Djakov für 


10) E. Djakov, Hochfrequenztechn. u. Elektroak. 58, 54 (1939). 
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parallele Leitwerte erhalten wir: 


_t-(1— 
b-(1 + #) 
(7) 


In diesen Gleichungen ist t = tg (a !), wobei / die Resonanzlänge der mit & 
abgeschlossenen Meßleitung darstellt. Weiterhin gilt: 


+ Stet 8 
4 tg? (« x) (8) 


Darin ist x die Länge, um die die Brücke K aus der Resonanzstellung ver- 
schoben werden muß, damit der Strom eines Indikators mit quadratischer 
Kennlinie auf die Hälfte des Resonanzwertes sinkt. Die Größe b wird durch 
die Beziehung 

. 
für den kleinsten Eingangswiderstand der MeBleitung, also für den Resonanz- 
widerstand, charakterisiert. In Abb. 3 09 


ist b in Abhängigkeit von + aufge- ! 


tragen. Bei kleinen I kann man in % 
(8) entwickeln und erhält 


b=ar für (10) 


Abb. 3. b als Funktion von = 
trägt. A 

Aus den von Djakov aufgestellten Beziehungen kann man weiterhin 
folgern, daß der Betrag des Eingangswiderstandes der MeBleitung bei einer 


Verstimmung der Brücke um rd gegen die Resonanzlage den Wert R (i + $) =Z 


hat. Mit (9) ergibt sich dann für einen Gleichrichter mit quadratischer Kenn- 
linie : 


_ Emm _ 


(11) 


Dieser Weg, b zu bestimmen, wird bei großer Dämpfung durch die Probe vor- 
zuziehen sein. 


Es wurde schon gezeigt, wie die Gleichungen für e’ und e” sich im Falle 


9 = 0 vereinfachen. Aus (6) folgt, daß p = 0 eintritt für ¢ = 0 oder/ = n 2 
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und für i= co oder] = (2 + „2 7 (n ganze Zahl). Unter Berücksichtig 


der Transformationseigenschaften der am Anfang kurzgeschlossenen Mei, 
leitung entnimmt man der Abb. 2, daß t = 0 der Länge s,, der Probenleitung 


t = co dagegen der Länge s, entspricht. Wir setzen nun | = (2n + 14 


voraus. Die zugehörige Länge der Probe wäre im Falle &, = 0 gleich eine 
halben Wellenlänge in der Flüssigkeit. Bilden wir in (6) und (7) den Lime 
für = oo, so erhalten wir: 


9=0%; g=b für I-@n+nt. 


Da der Eingangsleitwert der Probe für konstanten Abschlußleitwert nur ei 
Funktion von ihrer Länge ist, sind / und 6 gleichfalls allein von s abhängig 
Durch Differentiation von (6) und (7) folgt: 


dp 1+® 
ds + 
dg 
ds 


fa - (t 


Bildet man hier ebenfalls den Grenzwert fiir t = oo, so ergibt sich: 


für 1=(2n+1)4. 


Diese Beziehungen gelten exakt, sie können jedoch durch Vernachlässigung 
noch vereinfacht werden, wenn die Voraussetzungen für (10) erfüllt sind. Mi 
(x x)?< 1 können wir dann schreiben: 


ais 


(18) 


d 


Aus den Gleichungen können wir entnehmen, daß die Eigenverluste der Lei 
tung und die Kopplungsdämpfung einen relativ geringen Einfluß auf das Meß 
ergebnis haben, da eine konstante Zusatzdämpfung z.B. in die Formel 
für e’”’ nicht eingeht. 
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3. Beschreibung der Meßanordnung 


Der Meßkreis ist in Abb. 4 im Schnitt dargestellt. Der Innendurchmesser 
des starkwandigen, geschliffenen Messingrohres, das als Außenleiter dient, 
verhält sich zum äußeren Durchmesser des konzentrischen Innenleiters wie 
3,8:1. Der Innenleiter besteht aus einem a 
relativ dünnwandigen Rohr, das am unteren N X "ihrometertrieb 
Ende offen ist. Die beiden Leiter sind im N far den Innenieiler 
oberen Teil durch einen gut eingepaßten 
Kolben kurzgeschlossen, der innen und 
außen Federkränze besitzt. Kurzschluß- 
kolben und Innenleiter können mit Hilfe 
von zwei Mikrometertrieben unabhängig 
voneinander bewegt werden. Der Über- 
sichtlichkeit wegen wurden die Koppel- AN: Mikrometertrieb 
schleifen, die am Kurzschlußkolben be- Ne für Kurzschlußkolben 
festigt sind, in Abb. 4 nicht mit eingezeichnet. 
Die Ebenen der Schleifen stehen ungefähr 
senkrecht aufeinander, damit eine direkte 
Kopplung vom Generator auf den Gleich- 
richterkreis vermieden wird. Durch Variation 
der Eintauchtiefe kann der Kopplungsgrad 
verändert werden. Die Zuleitungen zu den 
Schleifen werden innerhalb des Meßkreises 
von Messingrohren gebildet, die einen kon- 
zentrischen Innenleiter enthalten und außen Kurzschlußkolben 
ineinem elektrisch angepaßten Verbindungs- \ 


stück enden. Diese Anordnung war not- 
wendig, weil bei einem unmittelbaren An- 
schluß von flexiblen Kabeln Störungen auf- i 7 
traten. 


+ Steigrahr 
Der untere Teil des Meßkreises ist von 
einem Thermostatengefäß umgeben und kann 
vom oberen Teil getrennt werden. Seine 
Länge ist so gewählt, daß bei der tiefsten 
Stellung des Innenleiters der Abstand des 
unteren Endes von der Abschlußplatte noch 
größer ist als der Leitungsdurchmesser. u 
Damit soll erreicht werden, daß für Wellen- 
längen, bei denen keine Hohlraumresonanzen 
> unteren Teil auftreten, der Abschluß- 
itwert der Probenleitung nicht von der , 
Stellung des Innenleiters abhängt. 4. 


Die Meßflüssigkeit tritt durch eine Bohrung in der Abschlußplatte in die 
Leitung; ihre erforderliche Höhe wird mit Hilfe einer feinen Marke im Steig- 
rohr eingestellt. Der ganze Meßkreis ist auf ein schweres, justierbares Stativ 
montiert, das die nötige Stabilität gewährleistet und eine senkrechte Einstel- 
lung der Leitung ermöglicht. Als Generator wurde ein Leistungsmeßsender 
mit der LD12 in Gitterbasisschaltung verwendet, dessen Netzteil aus einem 


Ann. Physik. 6. Folge, Bd. 17 4 
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Netzspannung-Konstanthalter gespeist wurde. Die Konstanz der Wellenläng 
und der Hochfrequenz-Amplitude wurde mit Hilfe eines Absorptionskreisg 
während der Messungen überwacht. 


Zwischen Generator und Meßkreis ist ein beidseitig angepaßter Ohmsche 
Spannungsteiler 1:6 bzw. 1:16 eingeschaltet, wodurch Rückwirkungen auf 
den Generator vermieden werden. Die Verbindungen sind durch flexibk 
Koaxialkabel vom Wellenwiderstand 70.2 hergestellt, die zur zusätzliche 
Abschirmung mit Kupferfolie bewickelt wurden. Hinter der Auskoppel 
schleife wird die Hochfrequenzspannung in einem Siliziumdetektor vom 
Typ ED 704 (OSW) gleichgerichtet. Der Richtstrom wird mit einem ebenfalk 
sorgfältig abgeschirmten Lichtmarkengalvanometer gemessen. 


4. Durchführung der Messungen 


Zunächst wird die Kennlinie J = const i? des Indikators ermittelt, indem 
am Topf ohne Probe eine Resonanzkurve aufgenommen und ihr Verlauf mit 
dem theoretischen verglichen wird. Dieser Vergleich zeigt auch, ob wirklich die 
vorausgesetzte Lecherschwingung allein angeregt wird. Djakov!P) leitete 
für das Quadrat des auf den Resonanzwert bezogenen HF-Stromes am Eingang 
der MeBleitung ab: 


(AD) + tg? (a - Al) 
ies 1 + tg? (a Al) 


Al ist dabei die Verstimmung gegen Resonanz. Für b?< 1 und a - Al<1 
folgt 


1— 1 


Trägt man nun : = als Funktion von Al im doppellogarithmischen Maß- 


stabe auf, so ergibt sich eine Gerade, falls man die richtige Kennlinienpotem 
des Gleichrichters angenommen hat und falls die Anregung einwandfrei ist. 
Gleichzeitig können wir hier den Fehler abschätzen, der durch Annahme eine 
falschen Kennlinie in 6 auftritt. Für y=} ist Al gleich der Halbwertver- 
stimmung x. Ist die fälschlich angenommene Kennlinienpotenz q, die wirk 
liche dagegen p, so mißt man x bei einem Stromwert, der y = (4)*/? entspricht. 
Bei kleinem x ist demnach der gemessene Wert von b um den Faktor y 2:!/? —1 
falsch. Nimmt man z. B. lineare Gleichrichtung an (q = 1), liegt aber qua 
dratische vor (p = 2), so wird b um 36% zu klein gemessen. Bei größerem # 
sind die prozentualen Fehler etwas kleiner. Die Überlegung zeigt also, dab 
beachtliche Fehler in b und damit vor allem in e” auftreten, wenn die Indikator 
kennlinie nicht genau bekannt ist. Deshalb ist es zweckmäßig, die Kennlinie 
auch noch dadurch zu bestimmen, daß man vor den Meßkreis nacheinandet 
Spannungsteiler von verschiedenen Teilungsverhältnissen einschaltet und 
damit die eingekoppelte Spannung in definierter Weise variiert. 


Ferner ist zu prüfen, ob der Abschlußleitwert der Probenleitung bei dem 
verwendeten Wellenlängen wirklich unabhängig von der. Stellung des Innen 
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leiters ist. Zu diesem Zweck wird der Meßkreis ohne Flüssigkeit in Resonanz 
bracht; sodann wird der Innenleiter um eine bestimmte Strecke verschoben 
und schließlich der Kreis mit Hilfe des Kurzschlußkolbens wieder abgestimmt. 
Ist die Streukapazität konstant geblieben, so ist die elektrisch wirksame 
Länge des Innenleiters bei konstanter Resonanzwellenlänge stets die gleiche. 
In der Flüssigkeit kann die Prüfung durchgeführt werden, indem man die 
Flüssigkeitshöhe um einen bestimmten Betrag ändert und den Innenleiter ent- 
sprechend nachstimmt. 


Die eigentliche Messung geht in folgender Weise vor sich. Zunächst wird 
die MeBleitung auf / = J eingestellt. Sie wird deshalb statt mit der Proben- 


leitung mit einer Kurzschlußscheibe abgeschlossen und so in eine + -Resonanz- 


leitung verwandelt. Wird nun mit dem Kurzschlußkolben auf Resonanz abge- 
stimmt, so ist die Meßwellenlänge ermittelt. Wird weiterhin die Meßleitung um 


dieHälfte ihrer Länge verkürzt, so ist die gewünschte 4. -Binstellung erreicht. 


Jetzt kann die MeBleitung wieder mit der Probenleitung abgeschlossen werden, 
in die soviel Flüssigkeit eingefüllt wird, daß deren Oberfläche gerade dort 
liegt, wo vorher die Kurzschlußebene lag. Das läßt sich mit Hilfe von Marken 
oder Überläufen leicht erreichen. Stimmt man nun mit dem Innenleiter auf 
Resonanz ab, so hat die Resonanzkurve auch bei größeren dielektrischen Ver- 
lusten in der Flüssigkeit noch relativ steile Flanken. In der Umgebung von 


I= + bzw. s= s, werden dann einige Wertetripel s, /, x bestimmt. Hier 
kann die Messung der geforderten Genauigkeit angepaßt werden ;in den meisten 
Fällen genügen zwei Meßpunkte neben sy. Von Vorteil ist, daß nach der 2 
Einstellung nur noch Längendifferenzen gemessen werden. 


Es ist unter Umständen zu berücksichtigen, daß die durch die Innenleiter- 
wand verdrängte Flüssigkeitsmenge sich mit s ändert. Ist der Querschnitt der 
Rohrwand nicht sehr klein im Verhältnis zum Querschnitt der Leitung, so 
muß an As die entsprechende Korrektur angebracht werden. 

Aus den Meßwerten ergeben sich die Differenzenquotienten = und = 
und die Größe b(s,,), mit deren Hilfe man nach den Gln. (16) und (17) e’ und e” 
errechnet. Falls bei geringer Dämpfung die Beziehungen (18) verwendet 


werden, genügen = und = . Die Genauigkeit der Ergebnisse hängt vor 


allem davon ab, wie genau die Differentialquotienten z und = bestimmt 


werden. Eine Abschätzung des Fehlers, der entsteht, wenn man die Diffe- 

renzenquotienten einsetzt, ist mit Hilfe der Krümmungen at und as 
an der Stelle s,g möglich. Wenn, wie es bei der beschriebenen Anordnung der 
Fall war, die Senderintensität auf 0,5% und die Meßfrequerz auf 0,1% kon- 
stant bleiben, ergeben sich aus der Abschätzung für A = 30cm und für s- 
Intervalle von As = 0,05cm maximal folgende prozentuale Fehler von e’ 
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5. Einige MeBergebnisse 


Nach der beschriebenen Methode wurden unter anderen Messungen an 
binären Flüssigkeitsgemischen (Alkohol und eine unpolare Komponente) und 
an ternären Systemen (zwei verschiedene Alkohole und eine unpolare Kom- 
ponente) durchgeführt. Aus der Abhängigkeit der dielektrischen Konstanten 
von der Temperatur und von der Konzentration sind Aussagen über die Asso- 
ziation bzw. Mischassoziation zu gewinnen. Die Ergebnisse sollen an anderer 


Tabelle 1 Tabelle 2 
Methylalkohol-Benzol Isobutylalkohol-Benzol 
Freoput = 0,1617 


#°C 


20 
30 
40 
50 
60 


Tabelle 3 
Methylalkohol- 
Benzo 


= 0,0769, = 0,0772 
°C 


Stelle ausführlich veröffentlicht werden. Hier sei nur als Beispiel eine Meß- 
reihe angeführt, die die Untersuchungsmethodik kennzeichnet. 

Als unpolare Flüssigkeit wurde Benzol gewählt, weil sich aus Dipolmoment- 
Bestimmungen bei 3MHz ergab, daß sich dieses Lösungsmittel gegenüber 


52 ee Th. 
und e”: Methy 
| * (%) im Sit 
| | dann, 
von de 
10-2 | 0,4 3,5 
100 | 10% | 05 4,8 
| | 8’ | 8 || | €’ | 
20 3,14 0,237 3,02 0,079 
30 2,97 0,138 2,91 0,049 
| 40 2,91 0,061 2,79 0,023 
50 2,80 0,029 932 | 0,012 
: 60 2,79 0,008 2,56 | 0,004 
20 2,94 0,137 : 
30 2,79 0.096 =~ 
| 40 2,77 0,054 i 
50 2,70 | 0,032 güns 
60 247 | 0,012 Lösu 
berei 
: den 
Benz 
1:1; 
gefül 
I 
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Methylalkohol und Isobutylalkohol, d. h. gegenüber den beiden verwendeten 
Dipolkomponenten, ungefähr ebenso milieugleich verhält. Ideal milieugleich 
im Sinne Euckens sind eine Dipolflüssigkeit und eine unpolare Flüssigkeit 
dann, wenn sich die zwischenmolekularen Kräfte zwischen den Dipolmolekülen 
von den Kräften zwischen Dipolmolekülen und unpolaren Molekülen nur durch 


Mettylalkohol-Benzol 


[(els] 60° 


Abb. 5. Realteil und Imaginärteil der Di- 

elektrizitätskonstanten von Alkohollösungen 

in Abhängigkeit von der Temperatur bei 
A = 30 cm 


die Dipolkräfte unterscheiden, die im ersteren Fall zusätzlich vorhanden sind. 
Diese Bemerkung ist hier notwendig, weil in bezug auf die Milieugleichheit 
bei Methylalkohol z. B. Tetrachlorkohlenstoff als unpolare Komponente 
günstiger wäre. In Tabelle 1 sind die Meßergebnisse bei A = 30 cm für eine 
Lösung von Methylalkohol (Molenbruch fyetn) in Benzol im Temperatur- 
bereich von 20° C bis 60° C verzeichnet. Tabelle 2 gibt die entsprechen- 
den Werte für eine Lösung von Isobutylalkohol (Molenbruch fisoput) in 
Benzol. Schließlich wurden die beiden Lösungen im Volumenverhältnis 
1:1 gemischt und die dielektrischen Messungen am ternären System durch- 
geführt (Tabelle 3). 


Die Temperaturangaben besitzen eine Genauigkeit von +0,5° C. 
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In Abb.5 ist zum Vergleich der Verlauf von e’ und e” bei den dei 
Lösungen dargestellt. 


Herrn Mechanikermeister H. Preußer möchte ich für die präzise Arbeit 
beim Bau der Meßleitung danken. Ferner danke ich Fräulein D. Dietze, 
die mir bei den Messungen hervorragende Hilfe leistete. 


Leipzig, Physikalisches Institut der Karl-Marx-Universität. 


Bei der Redaktion eingegangen am 12. April 1955. 
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Eine einfache Methode, magnetische Elementarbezirke 
mit trockenem Pulver sichtbar zu machen 


Von W. André und E. Schwabe 
Mit 2 Abbildungen 


Bei der bekannten Bitterstreifentechnik werden im allgemeinen zur Sicht- 
barmachung der magnetischen Elementarbezirke kolloidale Lésungen von 
magnetischen Teilchen in Flüssigkeiten verwendet. H. J. Williams, R. M. 
Bozorth und W. Shockley!) haben bereits 1949 auch mit trockenem 
Pulver gearbeitet. Sie lieBen ein Aerosol, das Karbonyleisenpulver mit sehr 
kleinen Teilchen (Durchmesser kleiner als 1 «) enthielt, auf einer elektroly- 
tisch polierten Eisen-Silizium-Probe absitzen und erhielten Pulvermuster, 
die mit den Mustern der flüssigen Suspension übereinstimmten. 

Da die Gewinnung sehr feinen Karbonyleisenpulvers Schwierigkeiten 
macht, wurde von uns die folgende sehr einfache Methode zur Herstellung 
eines geeigneten Aerosols angewendet: Unter dem Abzug wurde Eisenpenta- 


Abb.1. Magnetische Elementarbereiche auf Kobalt, 
mit trockenem Fulver sichtbar gemacht. V = 150 


karbonyl in eine flache Schale gegossen und angezündet. In den entstehenden 
Rauch, der sehr viele feine Teilchen eines magnetischen Eisenoxyds enthielt, 
wurde die in der üblichen Weise vorbereitete Probe (in diesem Falle Kobalt) 
gebracht. Die Abb. 1 und 2 zeigen die auf diese Weise gewonnenen und die 


) H. J. Williams, R.M. Bozorth u. W. Shockley, Physic. Rev. 75, 155 (1949). 


AS 


Annalen der Physik. 6. Folge. Band 17. 1955 


mit einer flüssigen Suspension erhaltenen Pulvermuster. Wie zu erkennen is, 
werden mit dem trockenen Pulver beträchtliche Feinheiten wiedergegeben 
wenn auch das Auflösungsvermögen der flüssigen Suspension nicht erreicht 
wird. Bei geeigneter Versuchsführung ist jedoch eine Verbesserung der Methoi 
durchaus zu erwarten. 


Abb. 2. Wie Abb.1. Zur ‘Sichtbarmachung un 
eine kolloidale Lösung von Fe,O, in Seifenwasser 
verwendet 


Die Methode des trockenen Pulvers hat den Nachteil, daß die erhaltene 
Muster unbeweglich sind. Andererseits bietet sich hierbei der Vorteil, bei 
tiefen und hohen Temperaturen arbeiten zu können. 


Jena, Institut für N Werkstoffe der Deutschen Akademie der 
Wissenschaften zu Berlin. 


Bei der Redaktion eingegangen am 12. April 1965. 


Verantwortlich 
für die Schriftleitung: Prof. Dr. Friedrich Möglich, Berlin-Buch, Lindenberger Weg 10; 
für den Anzeigenteil: VEB Georg Thieme, Anzeigenabteilung, Leipzig C 1, Hainstr. 17-19, Aufg.C, 
Ruf 21981. Z. Z. gilt Anzeigenpreisliste Nr. 4; Verlag: Johann Ambrosius Barth, Leipzig Ol, 
Salomonstr. 18B, Fernruf: 63105, 63781. Veröffentlicht unter der Lizenznummer 285/12 
des Amtes für Literatur und Verlagswesen der Deutschen Demokratischen Republik 
Printed in Germany Druck: Paul Dünnhaupt, Köthen (IV/5/1) L 1881 
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Praxis der Staubmessung 
Ein Leitfaden der Meßtechnik mit 


dem Zeiß-Konimeter 


Von KÄTE GÜNTHER, Jena 


V, 42 Seiten mit 26 Abbildungen im 
Text. 1954. DIN A5. DM 4.30 


JOHANN AMBROSIUS BARTH/VERLAG/LEIPZIG fiir Wissenschaft und Technik 


Zeitschrift 
fiir wissenschaftliche Photographie, 
Photophysik und Photochemie 


Unter Mitwirkung von 
Prof. Dr. H. Frıeser, Köln — Prof. Dr. E. STEnGER, Kreuzwertheim a. M. 
herausgegeben von Dr. H. Arens, Dessau 


Jubiläumsband 50 
Teil I und II (1955) 


Für diesen Jubiliumsband, der Ende November 1955 erscheinen wird, sind von 
Fachvertretern des In- und Auslandes so zahlreiche Arbeiten zur Verfügung ge- 
stellt worden, daß es nicht möglich ist, sie im Rahmen eines normalen Bandes 
unterzubringen. Deshalb werden beide für das Jahr 1955 vorgesehenen Bände 
unter der Bezeichnung Band 50 Teil I und II, von denen jeder einem normalen 
Bandumfang entspricht (Preis je DM 28.—), ausgegeben. 

Die Zeitschrift ist eine Sammelstätte für alle Arbeiten auf wissenschaftlich- 
photographischem Gebiete, die Aufklärung photographischer Vorgänge im wei- 
testen Sinne bezwecken. Sie wendet ihr Interesse allem zu, was mit der Photo- 
physik und Photochemie zusammenhängt. Ebenso kommen entsprechend der 
Bedeutung der Farbenphotographie für das moderne wissenschaftliche und 
kulturelle Leben auch Gebiete der wissenschaftlichen Farbenlehre zu Wort, so- 
weit sie mit der Farbenlehre eng verknüpft sind. 


JOHANN AMBROSIUS BARTH / VERLAG / LEIPZIG 


Praktische Finführung in die physikalische Chemie 


Von Prof. Dr. K. LOTHAR WOLF, Kirchheimbolanden 
und Dr. habil. HANS GEORG TRIESCHMANN, Ludwigshafen/Rh. 


3., neubearbeitete Aufla 
VII, 255 Seiten mit 66 Abb. i. Text. 1954. Gr. 8°. DM 8.10, Halbleinen DM 9.70 


Tonindustrie-Zeitung u. keram. Rundschau: Auf leicht verständliche und trotz- 
dem wissenschaftliche Art vermitteln die beiden. Verf. zuerst die theoretischen Grundlagen 
und daran anschließend die praktischen Versuche, die zum Verständnis der Physikalischen 
Chemie geeignet sind. Besonders erfreulich ist die klare, einfache Behandlung des Stoffes, 
zu besonders angenehm bei den schwierigeren elektrochemischen Problemen ins Auge 


JOHANN AMBROSIUS BARTH / VERLAG / LEIPZIG 
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Statistik und Physik der Kometen 


Von Dr. NIKOLAUS B. RICHTER 
Sternwarte der Deutschen Akademie der Wissenschaften zu Berlin in Sonneberg (Thür)# 


VI, 142 Seiten mit 58 Abbildungen im Text. 1954. Gr. 8°. DM 13.40 


Die Sterne: Als langjähriger Spezialist auf dem Gebiet der Kometenforschung versteht 
N. B. Richter es meisterhaft, in allen Fragen das Wesentliche hervorzuheben und den Kern 
der Probleme aufzuzeigen. Das Buch ist unter strenger Wahrung wissenschaftlicher Exakt 
heit so gemeinverständlich und zugleich fesselnd abgefaßt, daß es auch jedem Sternfreund™ 
zur Anschaffung empfohlen werden kann. Druck und Ausstattung des Werkes sind vor 
züglich. M. Beyer 


Meteorströme 
Meteorie Currents 


Von Prof. Dr. C. HOFFMEISTER 
Sternwarte der Deutschen Akademie der Wissenschaften zu Berlin in Sonneberg (Thür,) 


VIII, 286 Seiten mit 109 Abbildungen im Text und 1 Ausklapptafel. Mit Zu-@ 
sammenfassungen in englischer Sprache. 1948. Gr. 8°. DM 18.—, Hin. DM 20.—9 


„Orion‘: Das Buch wendet sich in erster Linie an die Fachleute, und zwar Astronomen, 
Meteorologen und Geophysiker, die daraus Aufschluß erhalten sowohl über rein astro- 
nomische Fragen, wie die Beziehung der Meteore zu Kometen, als auch über gewisse 
terrestrische Probleme, die mit dem Zufluß kosmischer Materie zu den oberen Luftschichten 
zusammenhängen. Zweitens wendet es sich an die Sternfreunde im allgemeinen, die sich 
diesem Zweig der Forschung, sei es praktisch oder nur zu Studienzwecken, widmen wollen, 
R.A.NZG 


Temperaturmessung 
Von Prof. Dr. F. HENNING, Bremen 


2., umgearbeitete Auflage, VI, 304 Seiten mit 80 Abbildungen im Text. 1956, : 
Gr. 8°. DM 16.80, Leinen DM 19.— ; 


„Die Technik“: Das vorliegende Buch darf als eine Enzyklopädie des gesamten Ge 
bietes der Temperaturmessung in Laboratorien angesprochen werden. Als solches ist es 
unentbehrlich für jeden Physiker und Ingenieur. Ein Standardwerk auf dem Gebiet der 
Temperaturmessung. 
„Archiv für die gesamte Wärmetechnik‘‘: Das Buch ist außerordentlich flüssig geschrieben 
und geht auch auf thermodynamische Grundfragen ein (es enthält hier mehr an Grund“ 
legung, als Titel und Kapitelüberschriften vermuten lassen). Das Werk stellt wohl die 
beste heutige Darstellung des Gesamtgebietes dar, die für jeden Ingenieur und Physiker 
unentbehrlich ist, der sich mit Temperaturmessung und verwandten Fragen zu befassen hat, % 


JOHANN AMBROSIUS BARTH / VERLAG / LEIPZIG 
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